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Abstract 


The aim of this review is to consider optical characteristics of terrestrial clouds. Both single and multiple light-scattering 
properties of water clouds are studied. The numerous results discussed can be used for solutions to both inverse and direct 
problems of cloud optics. An introduction to main microphysical characteristics of water and ice clouds, including the effective 
size, the liquid water content and path, the concentration of particles, is given. The refractive index of water and ice, which is of 
importance for optical waves propagation in clouds, is also considered. We consider extinction and absorption coefficients, 
phase functions, and asymmetry parameters of water and ice clouds and their relationships with the size, shape, and 
concentration of water droplets and ice crystals. Simple analytical results for absorption, reflection, and transmission 
characteristics of cloudy media are given. They are valid for optically thick clouds with optical thickness larger than 5. We have 
reviewed studies of polarization of solar light reflected from clouds. Image transfer problems have been considered as well. 
Modern methods of passive remote sensing of water and ice clouds are discussed. The review is finished by the discussion of 
optical properties of inhomogeneous clouds, which is a main topic of modern cloud optics research. 
© 2003 Elsevier B.V. All rights reserved. 


Keywords: Radiative transfer; Light scattering and propagation; Water and ice clouds 


serve as a blanket to protect the Earth’s surface against 
cooling at night. This is due to the fact that the 
maximum of the terrestrial emission is in the far infra- 


1. Introduction 


Water, ice, and mixed clouds are major regulators of 


solar fluxes in the terrestrial atmosphere (Kondratyev 
and Binenko, 1984; Liou, 1992, 2002). They reflect a 
large portion of incoming visible radiation back to 
outer space. The light energy absorbed by water drop- 
lets and ice crystals leads to heating of atmospheric 
layers. Another important role of cloudy media is to 
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red, where water droplets are highly absorbing. 
Clouds, mists, and fogs are very common. They 
reduce the visibility of objects in the atmosphere and 
limit capabilities of atmospheric vision, remote sens- 
ing, and detection systems (Zege et al., 1991). Thus, it 
is of a great importance to understand the peculiarities 
of light interaction with cloudy media, which can be 
considered as a huge “colloid”, composed of liquid 
and solid water particles, dispersed in the air. It should 
be pointed out that water droplets and air around them 
can be contaminated by various fine particles (e.g., 
soot and dust particles). These aerosol particles influ- 
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ence both light-scattering and absorption properties of 
cloudy media, making the studies of light propagation 
in cloudy media even more difficult (Menon et al., 
2002). The very existence of clouds is due to fine- 
aerosol particles, the so-called cloud nuclei (Twomey, 
1977), with typical radii between 5 and 200 nm 
(Mason, 1975). 

Studies of optical and microphysical properties of 
cloudy media have a long and fruitful history. The 
results of these investigations have been summarized in 
numerous books and papers (Kondratyev and Binenko, 
1984; Zege et al., 1991; Liou, 1992, 2002; Kokhanoy- 
sky, 2001a). The main properties of cloudy media are 
well understood now, but unsolved problems remain. 
The most important of these include the accounting for 
the three-dimensional shape of clouds (Evans, 1998; 
Macke et al., 1999; Schreier, 2001; Schreier and 
Macke, 2001), and their inhomogenity in horizontal 
and vertical directions (Cahalan et al., 1994, 2001; 
Platnick, 2000, 2001). The influence of aerosols inside 
cloudy media, which can be one of explanations of the 
anomalous absorption paradox (Danielson et al., 1969; 
Rozenberg et al., 1978; Stephens and Tsay, 1990), also 
should be clarified. 

Another current issue is the characterization of 
optical properties of ice clouds, which have extremely 
complex microstructure (Liou, 1992, 2002; Liou and 
Takano, 1994; Macke et al., 1996; Yang and Liou, 
1998; Kokhanovsky and Macke, 1999; Yang et al., 
2000, 2001) and appear almost with the same fre- 
quency as water clouds. One rarely finds two identical 
crystals in ice clouds because preferential shapes vary 
with temperature and pressure (Mason, 1975), which 
vary across cloudy media (Yang et al., 2001). It 
should be stressed that optical properties of a single 
particle and, as a result, optical properties of a cloud 
as a whole are highly influenced by particles’ shapes 
(Mishchenko et al., 1995, 1996, 1999, 2000). 

A complete review across all fields of cloud optics, 
which is in a state of a rapid transition, would be 
difficult within a single book, not to mention in a single 
article. Thus, I will concentrate mostly on analytical 
results, which can be used for cloud optical properties 
studies. The formulae presented can be applied for 
rapid estimations of light fluxes in cloudy atmosphere 
and will assist understanding the information content of 
transmitted and reflection functions of cloudy media in 
respect to microphysical characteristics of clouds. 


2. Microphysical characteristics 
2.1. Water clouds 


2.1.1. Particle-size distributions 

Water clouds consist of small liquid droplets that 
generally have the spherical shape although particles 
of other shapes can exist due to different external 
influences. For instance, the deformation of large 
particles due to the gravitation force is of importance 
for raining clouds with particles having radii 1 mm 
and larger (Macke and Grossklaus, 1998). The aver- 
age radius of droplets in nonraining water clouds is 
usually around 10 um and the approximation of 
spherical particles works quite well. Natural clouds 
with droplets of uniform size throughout never occur 
due to the variability of physical properties of air both 
in space and time domains (Twomey, 1977). Thus, 
one can consider a radius of a droplet, a, as a random 
value, which is characterized by the distribution 
function fla). This function is normalized by the 
following condition: 


f flasa=1. (2.1) 
0 


The integral 


F(a) = / Naa (2.2) 


gives the fraction of particles with radii between a, 
and az in a unit volume of a cloud. The distribution 
function f(a) can be represented as a histogram, 
graphically or in a tabular form (Ayvazian, 1991). 
However, it is the most common to use an analytical 
form of this function, involving only two or three 
parameters (Deirmendjian, 1969). Of course, this is a 
great simplification of real situations occurring in 
natural clouds, but most optical characteristics of a 
cloud as a whole practically do not depend on the fine 
structures of particle-size distributions (PSD) f(a) 
(Hansen and Travis, 1974). McGraw et al. (1998) 
found that the local optical properties of polydisper- 
sions can be modelled with high accuracy by just the 
first six moments of the particle-size distribution, and 
the use of certain combinations of moments can 
reduce the number of parameters even further. 
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In most cases, the function f(a) can be represented 
by the gamma distribution (Deirmendjian, 1969): 


f(a) = Nate, (2.3) 
where 
pir! 
“Fare = 
H 0 


is the normalization constant and I'(u+1) is the 
Gamma function. It follows that the first derivative 
f'(aọ)=0 and the second derivative f"(ag) <0 (see Eq. 
(2.3)). Thus, the function (Eq. (2.3)) has the maximum 
at a=ao. Note that Eq. (2.4) follows from Eq. (2.1) 
and the definition of the Gamma function: 


T(u) = | = xt Lede. (2.5) 
0 

One can see that the parameter u characterizes the 

width of the particle-size distribution f(a), being 

smaller for wider distributions. Moments 


(a") = f tsaia (2.6) 


0 


of the distribution (Eq. (2.3)) are calculated from the 
following simple analytical equation 


m _ (9\"T(u+n+ 1) 
Ge a 27 


Eq. (2.7) can be used to find the average volume of 
particles 


(Vv) = +f aa de, (2.8) 


the average surface area 

(x) =4n 1 a’f (a)da, (2.9) 
0 

and the average mass of droplets 

(W) = p(V), 


where p=1 g/cm? is the density of water. It follows 
that 


(2.10) 


T(u+4) 

(V) aren (2.11) 
T(ut+3) 

(2) Sur" (2.12) 


T(u+4) 
== 2.13 
(W) PTF D ( ) 
where 
4 3 
vo = “ » = 4na;, wo = pvo (2.14) 


are corresponding parameters for monodispersed par- 
ticles with radii aọ. One can obtain in the case of the 
most often employed cloud model with aọ=4 um and 
u=6 (Cloud C.1 model of Deirmendjian, 1969): 


Ha a) = im, 
where vo = 2.7X107'° m3, so = 2x10? m’, 
wo =2.7x107'° g. Although parameters (Eq. 
(2.14)) are small, the very large numbers of cloud 
droplets (typically, 100 particles in cm?) create 
important factors for atmospheric processes (see next 
section). 

Eq. (2.3) allows to characterize the cloud droplet 
distribution only by two parameters: do and uw. 
However, it should be remembered that neither is 
constant and can vary inside the body of a cloud. 
Thus, their values depend on the averaging scale, 
with large averaging scales producing more broad 
particle-size distributions with smaller values of u. 
The value of u=2 was found to be rather repre- 
sentative (Khrgian and Mazin, 1952) and is advised 
to be used in low-resolution cloud satellite retrieval 
algorithms. The parameter “=6 (Deirmendjian, 
1969), used in the derivation of Eq. (2.15), is 
typical only for small averaging scales (Fomin 
and Mazin, 1998). General features of the droplet 
spectra in water clouds were studied in a great 
detail by, e.g., Warner (1973). 

Parameters ag and u are defined in terms of the 
specific unimodal cloud droplet distribution (Eq. 
(2.3)). It is of an advantage to characterize cloud 
particle-size distributions by their moments. 
Moments can be retrieved from optical measure- 
ments without reference to specific distribution laws 
(McGraw et al., 1998). 

Hansen and Travis (1974) found that the effective 
radius 


(2.15) 


(2.16) 


aef = 
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is one of the most important parameters of any 
particle-size distribution. It is proportional to the 
average volume/surface ratio of droplets. The param- 
eter (Eq. (2.16)) can be also defined for nonspherical 
particles as we will see later. The coefficient of 
variance (CV) of the particle-size distribution 


A => (2.17) 


(2.18) 


is also of importance, especially for narrow droplet 
distributions. The value of s is called the standard 
deviation. The CV, which is equal to the ratio of the 
standard deviation to the mean radius (a), is often 
expressed in percent. 

It follows for the PSD (Eq. (2.3)): 


3 1 
as = a9(14+2), A= (2.19) 
u I+ u 
and, therefore, 
1 1-4? 
Te > 2.20 
=P 14247" (een) 


For instance, we obtain at =3: aef=2ao, A=0.5, 
s=(a)/2. The effective radius, aes, is always larger 
than the mode radius, do. 

Eq. (2.20) gives the meaning of the parameter u in 
the PSD (Eq. (2.3)). In situ measurements show that 
the value of ap often varies from 4 to 20 um (Mason, 
1975) and we[2,8] in most of cases. It should be 
pointed out that clouds with smaller droplets are not 
stable due to the coagulation and condensation pro- 
cesses (Twomey, 1977). Particles with a large radius 
cannot reside in atmosphere for a long time due to the 
gravitation force. Thus, several physical processes 
lead to the existence of the most frequent mode radius 
range. One can obtain from Eq. (2.19) and inequality 
2<p<8 that the value of A&[0.3,0.6]. Thus, it 
follows that the standard deviation of the radius of 
particles in water droplets is from 30% to 60% of an 
average radius in most cases. Smaller and larger 
values of A do occur, but values of A smaller than 
0.1 were never observed (Twomey, 1977). Larger 


values of A may indicate the presence of the second 
mode in the range of large particles (Ayvazian, 1991). 

Eq. (2.19) and the results for ag and u just reported 
lead to the effective radius, aef, of water droplets being 
in the range from 5 to 50 um, depending on the cloud 
type. Near-global survey of the value of aef using 
satellite data shows that typically Sum < aef < 15 um 
(Han et al., 1994). We see that water clouds with 
der>15 um are rare. This can be used to discriminate 
satellite pixels with ice crystals even at wavelengths 
where ice and water absorption coefficients are almost 
equal. Such a possibility of discrimination is due to 
much larger (e.g., in 5—10 times) effective sizes of ice 
crystals as compared to droplets. The large size of ice 
crystals will reduce the reflection function in near 
infrared considerably as compared to droplets. This 
reduction can be easily detected. Note that clouds with 
aer>15 um are often raining (Masunaga et al., 2002). 
Pinsky and Khain (2002) showed that the threshold of 
the occurrence of drizzle is around aef= 15 um. Then 
strong collisions of droplets start. 

Some authors prefer to use the representation of the 
particle-size distribution by the following analytical 
form (Ayvazian, 1991): 


In? A 
1 a 
= m 221 
fla) V2n0a exp 202 ? ( ) 


which is called the lognormal distribution. The rela- 
tions between the values of des, (a), A and the 
parameters of the gamma and the lognormal par- 
ticle-size distributions are presented in Table 1. The 
value of Aef in this table represents the effective 
variance defined as 


f ET 
a% f "arida 


This value is often used instead of the coefficient of 
variance, A, (Hansen and Travis, 1974) because of a 
special importance attached to the value of the effec- 
tive radius of droplets, aes, in a cloudy medium as 
compared to the average radius, (a), of droplets. For 
instance, light extinction in clouds is governed mostly 
by values of aef and liquid water content independ- 
ently on the particle-size distribution (Kokhanovsky, 
2001a). 


Aef = 
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Table 1 
Droplet-size distributions 
fa) B (a) Get A Aer (a") 
sees -ut pitt ay\"T(1+ +n) 
Gamma distribution Bate “o ao(1 + (1/p) ao(1 + (3/4) l/u+1 1/⁄(u+3 
rad) o( W) a w) (1/u+1) 1u+3) u) TOFA 
B pa I 2 2 2 ha 
Lognormal distribution a exp ar ian ame?” Aer? ev —1 e — 1 ate? 


2.1.2. The concentration of droplets 

The number concentration of particles, N, in addi- 
tion to size and shape of particles, is of importance for 
the optical waves propagation, scattering, and extinc- 
tion in cloudy media. The concentration of droplets 
depends on the concentration Cy of atmospheric 
condensation nuclei. The value of Cy is smaller over 
oceans than over continents; thus, the concentration of 
droplets in marine clouds is on average smaller than 
over continents. Generally, the smaller concentration 
of droplets over oceans means that they can grow 
larger, producing clouds with larger droplets over 
oceans, which is confirmed by the analysis of the 
satellite optical imagery as well (Han et al., 1994). 
This influences the occurrence and the rate of precip- 
itation. 

Swensmark and Friis-Christensen (1997) and 
Marsh and Swensmark (2000) have speculated that 
cosmic ray ionization could influence the production 
of condensation nucleii and, therefore, cloud proper- 
ties. This is of importance for climate change prob- 
lems as well (Swensmark, 1998). 

The dimensionless volumetric concentration of 
droplets C,=N(V) and the liquid water content 
(LWC) Cy=pC, or (see Eq. (2.10)) Cw=N(W) are 
often used in cloud studies as well. 

The typical variability of N, C,, and Cw in water 
clouds is presented in Table 2 (Mason, 1975; Fomin 
and Mazin, 1998). Numbers in Table 2 are represen- 
tative, but it should be remembered that these values 
can change in broader range in real situations. The 
liquid water content is not constant throughout a 
cloud, but has larger values near the top of a cloud 
(Feigelson, 1981; Paul, 2000; Yum and Hudson, 
2001). 

The liquid water path (LWP), w, is defined as 


w= i C,,(z)dz, (2.22) 


where /=z.—z, is the geometrical thickness of a 
cloud. If C,,=const, 
w = Cyl. (2.23) 
The geometrical thickness of clouds varies, depending 
on the cloud type (Landolt-Bernstein, 1988). It usually 
is in the range 500—1000 m for stratocumulus clouds. 
Near-global data obtained by Han et al. (1994) from 
satellite measurements show that the liquid water path 
is typically in the range 50—150 g/m°. The annual 
mean is equal to 86 g/m” (Han et al., 1994). 

Cloud systems can easily cover area S = 10° km? 
(Kondratyev and Binenko, 1984); thus, the total 
amount of water W= wS (for idealized clouds having 
the same liquid water path for the whole cloudy area) 
stored in such a water cloud system is equal approx- 
imately to 10° kg, if we assume that w= 100 g/m’, 
which is a typical value (Han et al., 1994). This under- 
lines the importance of clouds both for climate prob- 
lems and human activity (e.g., crops production, etc.). 


2.2. Ice clouds 


Microphysical properties of ice clouds cannot be 
characterized by a single particle-size distribution 
curve as is the case for liquid clouds even if one 
considers relatively small volumes of a cloudy 
medium. This is due to extremely complex shapes 
of ice particles in crystalline clouds. Major shapes of 
ice crystals are plates, columns, needles, dendrites, 
stellars, and bullets (Mason, 1975; Volkovitsky et al., 
1984; Takano and Liou, 1989, 1995; Liou, 1992, 
2002; Yang et al., 2001). Combinations of bullets 


Table 2 

Typical range of values N, Cy, and Cw in water clouds 

N (om) G Cw (g/m) 
20-1000 1077-1076 0.01-1 
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and needles are also common. The Magano—Lee 
classification of natural crystals includes 80 shapes 
(Magano and Lee, 1966), ranging from the elementary 
needle (the classification index Nla) to the irregular 
germ (the classification index G6). 

The concentration of crystals, N, varies with 
height, often in the range 50—50,000 crystals per 
cubic meter. The ice water content 


G = NW), (2.24) 


where (W) is average mass of crystals, usually in the 
range 107—107 ' g/m?. Thus, the average mass of a 
crystal is in the range 2 x 10° 7-2 x 107° g. Crystals 
have a bulk density, p, less than that of the bulk ice 
(p =0.3—0.9 g/cm?) due to the presence of impurities 
and bubbles inside ice particles (Landolt-Bernstein, 
1988). The size of crystals is usually characterized by 
their maximal dimension, D, which is related to the 
effective size (Yang et al., 2000). D is usually in the 
range 0.1—6 mm for single crystals and 1—15 mm for 
snow crystal aggregates. Smaller crystals (e.g., with 
maximal sizes around 20 um) also are present in ice 
clouds (Yang et al., 2001). The mode Do of size 
distribution curves f(D) depends on the shape of 
crystals, with characteristic values of Do being 0.5 
mm for plates and columns, 1 mm for needles, 
sheaths, and stellars, and 2 mm for dendrites. Distri- 
bution curves can be modelled by gamma distribu- 
tions with half-widths of distributions, 41/2, being 
larger for larger values of D. Frequently, 41⁄2 = D, 
where D is the average size of crystals. 

Simple shapes of ice crystals (e.g., hexagonal 
prisms) can be characterized by two dimensions: the 
length of the prism, L, and the diameter, D = av3, 
where a is the side of a hexagonal cross section. Even 
in this most simple case, two-dimensional distribution 
functions, AD,L), should be used. Note that AD,L) can 
be approximately reduced to one-dimensional func- 
tions, AD), due to the existence of empirical relation- 
ships between the length of crystals and their diameter 
in natural clouds (Auer and Veal, 1970). For instance, 
for hexagonal columns with the semiwidth, ¢, and the 
length, L, (Mitchell and Arnott, 1994), ¢=0.35L at 
L<100 um and ¢=3.48,\/Z otherwise. Pruppacher and 
Klett (1978) give the following empirical relationship 
for plates: L=2.44830°4”* at C=[5 um, 1500 pm]. 
Similar relationships can be obtained for other shapes 
(Mitchell and Arnott, 1994). 


It should be pointed out, however, that the whole 
concept of a single-distribution function breaks down 
for ice clouds. Recall that the microphysical properties 
of most water clouds (as far as their optical properties 
are concerned) can be characterized by just three 
numbers (e.g., the effective radius of droplets, the 
coefficient of variance of the droplet distribution 
function f(a), and the concentration of droplets) (Deir- 
mendjian, 1969). 

By contrast, at least 80 multidimensional particle- 
size distributions are needed if one would like to use 
the classification of crystals developed by Magano 
and Lee (1966). This calls for the introduction of a 
new way for particle characterization in the case of 
complex particulate systems such as ice clouds. The 
same problem arises in the optics of mineral aerosol, 
blown from the Earth’s surface to atmosphere (Volten 
et al., 2001; Mishchenko et al., 2002). 

One of the possible solutions of the problem lies in 
the characterization of ice crystals in an elementary 
volume of a cloudy medium by the function 


N 


(a,b) =) eta) +) hE), 


r=1 


(2.25) 


where f(a) is the size distribution of particles of 
a regular shape (e.g., hexagonal plates or columns), 
fi(b) is the statistical distribution of particles with 
random surfaces or the so-called irregularly shaped 
particles. Note that instead of single variables (e.g., 
the radius of particles a), we need to introduce vector 
parameters g and b in the case of crystalline clouds. 
In particular, we have for the components of the two- 
dimensional vector parameter, a, with coordinates a, 
az in the case of an idealized cloud with hexagonal 
columns: a; =D, a.=L. 

Functions f(b) present statistical distributions of 
some statistical characteristics of particles (e.g., aver- 
age radii, correlation lengths, etc.). Values of c; and c,. 
give concentrations of different crystal habits. Clearly, 
the simplest case is to consider the function f(a, b) as 
a sum of two functions: 


Sf (@, B) = fila) + ofh), 


where function f(a) represents particles of a regular 
shape (say, hexagonal cylinders) and the function 
f(b) represents statistical parameters of a single par- 
ticle of irregular shape. 


(2.26) 
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This irregularly shaped particle can be presented, 
e.g., aS a fractal (Macke and Tzschichholz, 1992; 
Macke et al., 1996; Kokhanovsky, 2003a). It should 
be pointed out that the function (5) in this case 
represents “‘fictive’’ particles (Kokhanovsky, 2001a) 
which do not exist in a cloud at all. However, ice 
cloud optical characteristics, calculated using f 
(b), indeed represent quite well the optical character- 
istics of particles with extremely diverse shapes 
(Macke et al., 1996; Kokhanovsky, 2003a). 

A similar approach to optical characterization of 
irregularly shaped particles was developed by Pelto- 
niemi et al. (1989), Peltoniemi (1993), and Muinonen 
et al. (1996). It is based on Monte-Carlo calculations 
of light scattering by a large particle with a rough 
surface. The model of spheres with rough surfaces 
was successfully applied to the optical characteriza- 
tion of irregularly shaped aerosol particles (Volten et 
al., 2001). This model presumably can be extended for 
the case of ice clouds as well. For this, the parameters 
of the irregularity of the so-called fictive particle 
(Kokhanovsky, 2001a) should be changed accord- 
ingly. This is due to different morphology of ice 
crystals (Magano and Lee, 1966) as compared to 
mineral aerosol (Okada et al., 2001). Also note that 


real part of the refractive index 


the refractive index of aerosol in visible is larger than 
that for ice. 


2.3. Refractive index of liquid water and ice 


The complex refractive index of particles suspended 
in the atmosphere is another important parameter in 
atmospheric optics studies (Liou, 1992). This is due to 
the fact that not only the size, shape, and concentration 
of particles influence the light propagation in atmos- 
phere. The chemical composition, thermodynamic 
phase, and temperature of particles are of importance 
as well. The refractive index of water droplets and ice 
crystals varies with the temperature and has been 
tabulated by many authors (see, e.g., Hale and Querry, 
1973 for liquid water and Warren, 1984 for ice). 

Spectral dependencies of real and imaginary parts 
of the complex refractive index of water and ice are 
presented in Figs. | and 2. The differences in light 
absorption by liquid and solid water are considerable, 
a fact that can be used for the retrieval of the cloud 
thermodynamic phase (liquid, ice, or mixed-phase 
clouds) from ground (Dvoryashin, 2002) and satellite 
(Knap et al., 2002) measurements. In particular, we 
see that the absorption band around the wavelength, 


wavelength, um 


Fig. 1. Real part of the refractive index of water and ice. 
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2=1.55 um, for ice is shifted to larger wavelengths as 
compared to liquid water. This shift in the absorption 
band position can be easily detected with the use of 
modern spectrometers (see, e.g., Bovensmann et al., 
1999). 

The real part of the refractive indices of water, ny, 
and ice, n;, do not vary considerably in the visible and 
near-infrared regions of electromagnetic spectra. Gen- 
erally, n(A)<ny(/), with ny in the range 1.33-1.34 
for 2=0.4—1 um and n; in the range 1.30—1.32 within 
the same spectral band. Somewhat larger values of 
refractive indices occur at shorter wavelengths. 

The spectral variability of the imaginary part of the 
refractive index of water, which is responsible for the 
level of absorption of solar radiation by clouds, is 
much higher (see Fig. 2). It changes six orders of 
magnitude in the spectral range 0.4—2 um both for 
liquid water and ice. 

Different impurities in water droplets (mainly soot 
(Markel, 2002) and various aerosol particles (Twomey, 
1977)) can change the imaginary part of the refractive 
index of droplets (especially in visible, where water is 
almost transparent). This may influence the accuracy of 
modern cloud remote-sensing techniques (Nakajima et 
al., 1991; Zege et al., 1998c; Schuller et al., 2000; 
Asano et al., 2001). Note that Asano et al. (2001) found 
that the temperature dependence (Kuo et al., 1993) of 


liquid water absorption can also influence the cloud 
droplets’ size retrievals. 


3. Local optical characteristics 
3.1. Water clouds 


3.1.1. Extinction coefficient 

The information presented in Section 2 can be used 
as an input for studies of light interaction with cloudy 
media on a global or local scale. In particular, the 
attenuation of a direct light beam with the intensity, Jo, 
in a cloudy medium is governed by the following 
equation: 


I = h exp(—t/cosvo), (3.1) 


where Jp is the intensity of an incident light, Z is the 
intensity of the transmitted direct light, and Wo is the 
solar zenith angle, 


ge l eee (3.2) 


is the optical distance in a cloud, H is the geometrical 
distance, and Gext is the extinction coefficient. It is 
assumed that the cloud can be represented by a 


1 E 4 T i T 
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Fig. 2. Imaginary part of the refractive index of water and ice. 
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horizontally infinite homogeneous plane-parallel 
layer, which is, of course, a great simplification 
of a real cloud field (Cahalan et al., 1994, 2001). 
The value -1/o,,, is called the photon-free path 
length and gives the average distance between 
photon-scattering events in a cloud. Knowing the 
average number of scattering events, ñ, it is easy to 
find the average distance, L=ñl, which photon 
travels in the medium before it escapes out of the 
scattering layer. We have for the average total time, 
which photon spends in a cloud: t=L/c, where c is 
the group light speed. 

The value of / is often in the range 10—200 m for 
water clouds; thus, the value of Gext is in the range 
0.005—0.1 m~', depending on the cloud type, and 
perhaps even larger for dense fogs and mists. The 
extinction coefficient, o, x, determines the meteorolog- 
ical range of visibility (Zege et al., 1991; Liou, 1992), 
which is defined as Sm=3.91/6ext. In particular, we 
have at oex,=0.1 m™!: Sm=39.1 m. In a cloudless 
atmosphere, this number is approximately 1000 times 
larger. 

A more advanced understanding of the image 
reduction by cloudy media can be achieved on the base 
of the image transfer theory (Zege et al., 1991; Zege 
and Kokhanovsky, 1994, 1995; Barun, 1995, 2000; 
Zuev et al., 1997; Katsev et al., 1998). This theory 
treats a cloud as a high-frequency spatial filter. The 
cloud greatly reduces high spatial frequencies of 
observed objects, causing fine features of objects to 
be lost due to multiple light scattering and photon 
direction randomization process. Clouds with larger 
particles have “better” optical transfer functions 
(higher values) as compared to clouds with very fine 
droplets, a fact that can be used for cloud microphysical 
parameters monitoring purposes (Zege and Kokhanov- 
sky, 1992; Zuev et al., 1997). 

The extinction coefficient varies with height in a 
cloud (see Eq. (3.2)) and, for a given value of the 
vertical coordinate, z, is given by 


Oext = N | fa)Coxda, (3.3) 


where N is the number concentration of particles. The 
value of the extinction cross section, Cext, is obtained 
from Mie theory valid for spherical droplets (Shifrin, 
1951; Van de Hulst, 1957; Kerker, 1969; Bohren and 


Huffman, 1983; Kokhanovsky, 2003b). It is known 
that (Shifrin, 1951) 


Ca = 5/2 (3.4) 


in a good approximation for water droplets in the 
visible range of the electromagnetic spectrum. Here X 
is the surface area of a droplet. From Eqs. (3.3) and 
(3.4)~ 


N(2) 
T 


Thus, the extinction coefficient depends on the prod- 
uct of the number concentration of particles and their 
average surface area. By use of 

Cy 

(V) 
where (V) is the average volume of particles (see Eq. 
(2.8)) and C, is the volumetric concentration of 
particles, we get from Eq. (3.5) 


BG, 


(3.5) 


Oext = 


ext — ; 3.7 
Oext Daet ( ) 
where we introduced the effective radius 
3(V) 
da =. (3.8) 
(2) 


The extinction coefficient decreases with the size of 
particles at C, = const. 

For constant Cy, the value of ext can be expressed 
via the liquid water content, Cw, which is often 
measured in cloudy media: 

3Co 
© 2paer’ 
where p=1 g/cm? is the density of water. It follows 
for typical values age=6 um and C,,=0.4 g/m? that 
Jext=0.1 m~}. This is rather dense cloud with the 
ratio Jp/[=e at the distance 10 m for light incident 
perpendicular to a cloud layer (see Eq. (3.1)). 

We should underline two peculiarities of water 
clouds, which follow from Eq. (3.9). First of all, the 
extinction does not depend on the wavelength, 2, and, 
second, it depends only on the ratio of liquid water 
content to the effective radius, ae. The first point leads 
to the conclusion that clouds do not change the spectral 
composition of a direct light beam in visible. Thus, 
spectral transmittance methods of particulate media 


(3.9) 


Oext 
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microstructure determination (Shifrin and Tonna, 
1993) cannot be applied in cloud optics. Another im- 
portant property, which follows from Eq. (3.9), is inde- 
pendence of the value of Sext on the type of particle-size 
distribution. Even the width of the PSD is of no im- 
portance at a given value of the effective radius of 
particles. This is the reason why the satellite methods of 
cloud microstructure determination (Arking and 
Childs, 1985; Nakajima and King, 1990; Han et al., 
1994) are concerned mostly with the retrievals of aef 
values and the liquid water path (see Eqs. (2.23) and 
(3.9) 


2 
Ww = z pat 


(3.10) 


and not the droplet-size distribution itself. 

The type of particle-size distribution in clouds is 
not possible to find from extinction measurements in 
the visible. On the other hand, there is an advantage in 
that one can calculate the extinction coefficient of 
cloudy media even with limited information on stat- 
istical properties of droplets and crystals. This is 
especially important for ice clouds as we will see later. 

Eq. (3.9) is valid only in the visible range of the 
electromagnetic spectrum and for comparatively large 
droplets. Calculations for longer wavelengths should 


— Mie theory 
approximation 


be performed with Mie theory although a simple cor- 
rection to Eq. (3.9) allows for its use in the near-infrared 
region of the electromagnetic spectrum (Kokhanovsky 
and Zege, 1995, 1997a,b; Kokhanovsky, 2001a): 


a 3Cy i+ v 
ext 2 pder (kas)? í 


where k=27// and v is the parameter which only 
slightly depends on the dispersion of the particle-size 
distribution. It follows that for most typical widths of 
particle-size distributions in clouds v = 1.1 (Kokhanov- 
sky and Zege, 1997b). We get from Eq. (3.11) that the 
extinction coefficient of cloudy media slightly in- 
creases with the wavelength in the visible and near 
infrared. The spectral dependence of Sext for the 
gamma PSD (Eq. (2.3)) at ae¢=4 um, u=6 is presen- 
ted in Fig. 3 for the spectral range 0.2-5 um. The 
data were obtained with simple Eq. (3.11) and Mie 
theory, assuming that Cy=0.1 g/m°?. The value of v 
was equal to 1.1. The accuracy of Eq. (3.11) 
decreases with wavelength/size ratio, but the error is 
smaller than 5% at 2<2 um and aef> 4 um. Smaller 
values of dep are of a rare occurrence in natural clouds 
(Han et al., 1994). Eq. (3.9) gives a constant value of 
0.375 m` ' in the case, presented in Fig. 3. We see 


(3.11) 


Aum 


Fig. 3. Extinction coefficient of a cloudy medium with water droplets, characterized by the PSD (Eq. (2.3)) at aee=4 um, w= 6. The value of Cw 


is 0.1 g/m’. 
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that Eq. (3.11) has a high accuracy and gives the 
correct spectral dependence the extinction coefficient 
(at least up to the wavelength equal to 1.5 um). 

The accuracy of Eq. (3.11) at larger wavelengths 
can be increased if one accounts for the dependence of 
the extinction cross section on the refractive index of 
particles, which could be done, for instance, in the 
framework of the Van de Hulst approximation (Acker- 
man and Stephens, 1987; Kokhanovsky, 2001a). 
Parameterizations developed by Mitchell (2000) and 
Harrington and Olsson (2001) can be used even out- 
side the geometrical optics limit. 

Clearly, the accuracy of Eq. (3.11) increases with 
aef. Thus, the case presented in Fig. 3 provides the 
maximal error of the approximation (3.11) for terres- 
trial clouds having aef>4 um. Typically, errors are 
even smaller. 


3.1.2. Absorption coefficient 

Clouds both scatter and absorb incident radiation. 
The probability of photon absorption, f, is defined by 
the following equation: 


Tabs 
p= at, (3.12) 
where Gabs is the absorption coefficient and Gex is the 
extinction coefficient (see, e.g., Eq. (3.11)). The value 
of f is close to zero in the visible and near infrared, 
which allows to obtain simple equations for reflection 
and transmission functions of cloudy media in terms 
of series 877, where j is an integer number (Rozen- 
berg, 1962, 1967; Van de Hulst, 1980; Minin, 1988; 
Kokhanovsky, 2002b). 
The absorption coefficient is defined as: 


Oabs = vf Caps f (a) da 
0 


where Cabs is the absorption cross section (Shifrin, 
1951). It follows for the scattering coefficient, Gsca = 
Oext — Fabs. The absorption cross section can be cal- 
culated from Mie theory for spherical particles. 
Another way to calculate the absorption cross section, 
which is valid for nonspherical particles as well, is to 
use the exact integral (Kokhanovsky, 2001a; Markel, 
2002): 


Cabs = = = | ney F)E* (FP) br, 
jel 


(3.13) 


where Eo is the electric field of an incident wave, Eis 
the electric field inside a particle, k=27/2 is the wave 
number, m=n — iy is the refractive index, V is the 
volume of a droplet or crystal. The first coarse 
approximation is to assume that the value of 
E@)E F) ~B(n) Bor )E0'() (3.15) 
due to the weakness of absorption. Thus, it follows in 
this case 


Caps = B(n)aV, (3.16) 


where «=4zy/2 and we assumed that a particle is 
uniform (n and y do not depend on the radius vector 7’). 
Clearly, we have that B— 1 as n— 1 for arbitrarily 
shaped particles (Kokhanovsky, 2001a). It follows 
from Eqs. (3.13) and (3.16): 


Cabs = B(n)aN(V) (3.17) 
or 
Cabs = B(n)aCy. (3.18) 


We have calculated the value B=o,,,/aC\, for the PSD 
(Eq. (2.3)) with w=6 and a.¢=4 um with Mie theory 
and presented it in Fig. 4 as the function of the wave- 
length. One can see that the value of B almost does not 
depend on the wavelength. This means that also there is 
only weak dependence of B on n(A), x(4), and Xep= kaef 
at 2 <1.8 um. For larger values of 4, Eq. (3.15) 
becomes invalid due to higher absorption of light by 
droplets or crystals. We have from Fig. 4: B = 5/3. 

The accuracy of approximation (3.18) at B=5/3 is 
better than 5% below 2=1.8 um for values of der in 
the range 4—6 um. For larger particles, the accuracy 
reduces due to the neccesity to account for terms 
proportional to g? and higher. Then Eq. (3.18) should 
be modified (see, e.g., Kokhanovsky and Macke, 
1997; Kokhanovsky, 2001a). 

Eq. (3.18) can be written as 


B(n)aCy 


Cabs = ——_—; (3.19) 
where Cw is the liquid water content and p is the density 
of water. One can see that the spectral variation of the 
absorption coefficient of liquid clouds approximately 
coincides with the spectral dependence of the absorp- 
tion coefficient of liquid water itself (see Fig. 2) at least 
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Fig. 4. Spectral dependence B(4), obtained for the same conditions as in Fig. 3. 


up to 1.8 um. The same applies to ice clouds, namely, 
Cabs ~ 42x(A)/A, where 7;(/) is the imaginary part of 
the refractive index of ice. This feature can be used for 
the cloud thermodynamic phase determination by 
remote-sensing techniques (see Fig. 2). Note that B in 
Eq. (3.19) depends on the shape of particles (Kokha- 
novsky and Macke, 1997). 

It should be pointed out that water droplets can 
collect absorbing particles from surrounding air. Liquid 
aerosol particles penetrate into droplets. They dissolve 
and change the value of « in Eq. (3.18). Soot particles 
can form a discontinuous layer on the surface of a 
droplet. All these effects will produce an enhancement 
of the absorption coefficient as compared to simple Eq. 
(3.18) (Prishivalko et al., 1984; Kokhanovsky, 1989; 
Asano et al., 2001). Note that the effects of aerosol 
particles (e.g., soot) can be modelled by inserting into 
Eq. (3.18) the value of 


N 
& = by + X Citi, 
i=1 


where ,, and g; are absorption coefficients of water and 
impurities, respectively, and c; is the relative concen- 
tration of impurities. The problem of black carbon 
influence on the light absorption in cloud was studied 
in detail by Chylek et al. (1996) and Liu et al. (2002). 


(3.20) 


The parameter of considerable importance in the 
radiative transfer problems in clouds is not the absorp- 
tion coefficient itself, but the ratio of the absorption 
and extinction coefficients. It follows from Eqs. (3.11) 
and (3.19) for this ratio 


B = Boas, (3.21) 


where =(2/3)B(n)(1 + vx >)" is close to 1. The 
value of f is obtained from the measurement of the 
reflection function of clouds in the infrared (Naka- 
jima and King, 1990). Thus, Eq. (3.21) can be 
used to find the effective radius of droplets 
(dep ~ P/a), which corresponds to a given value 
of p. The uncertainty in the value of « (see Eq. 
(3.20)) can introduce additional errors in the 
retrieved values of the effective radius of droplets 
(Asano et al., 2001). 

Let us estimate this error, neglecting the term 


vxor ” for the sake of simplicity. Then it follows: 


36 


te= Spi 


(3.22) 
and, therefore, (Ôaep'aef) =( — da/x). Aerosol particles 
trapped in droplets make the value of « larger. Thus, 
the effective radius, dg, of polluted clouds (as 
obtained from Eq. (3.22) without account for pollu- 
tants) appears to be smaller as compared to values 
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obtained from in situ measurements (e.g., using laser 
diffractometers (Kokhanovsky, 2001a)). This was 
found also experimentally (Nakajima et al., 1991). 
To correct this effect, one should use the value of « 
(Eq. (3.20)) in Eq. (3.22). The problem, however, is 
complicated by the fact that values of absorption 
coefficients and concentrations of impurities are not 
known a priori. On the other hand, one can formu- 
late the inverse problem in such a way that they can 
be derived from measured spectral reflectances 
simultaneously with the value of the effective radius 
of droplets (Zege et al., 1998c). 

Eq. (3.18) is limited to the case of weak absorption. 
It can be generalized to account for larger absorption 
of light by water droplets. We found, parameterizing 
Mie theory results, that 
Cabs = BaCy[1 — gaet]. (3.23) 
Eq. (3.23) transforms to Eq. (3.18) for small values of 
absorption («aef — 0). 

In conclusion, we present the spectral dependence 
of the probability of photon absorption in Fig. 5. Data 
were obtained using Mie theory and approximate 
solution, given by Eq. (3.23) with account for Eqs. 


(3.11) and (3.12). The effective radii of particles were 
equal to 4 and 16 um. The accuracy of the approx- 
imation is better than 10% at wavelengths smaller 
than 2 um (see Fig. 6). They are typically smaller than 
5% at å < 1.8 um. 


3.1.3. Phase function 

Until now, we have considered only the extinction 
and the absorption characteristics of cloudy media. 
However, clouds not only attenuate propagated sig- 
nals and absorb light but also scatter incident light in 
all directions. Generally, the probability of photon 
scattering by a droplet depends on the size of the 
droplet and the energy of the incident photon (or the 
light wavelength). 

Now we should characterize the probability of 
photon scattering in a given direction, specified by 
the scattering angle, 0 (0=0 corresponds to the for- 
ward scattering). It is known that the intensity of light 
scattered by water droplets is much larger in the 
forward hemisphere than it is in the backward direc- 
tions (Deirmendjian, 1969; Kokhanovsky et al., 
1998a). 

Let us introduce the probability dP of light 
scattering in the direction, specified by the vector 
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Fig. 5. Probability of photon absorption, obtained for the same conditions as in Fig. 3. Results for the effective radius, a.¢=16 um are also 


shown. 
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Fig. 6. Error of the geometrical optics approximation for the probability of photon absorption, obtained from data presented in Fig. 5 for 
effective radii 4 and 16 um. Data for the effective radius 6 um are also shown. 


Q inside the solid angle dQ. This probability will be 
proportional to the value of dQ/4z. Namely, we 
have: 


dQ 


dP = x(Q) 7, 


(3.24) 
where x(Q) is the coefficient of proportionality. It 
follows from Eq. (3.24): 


| x(Q) a =P 
4r An 


The value of P represents the probability of photon 
survival in the scattering process. It is equal to 1 if 
there is no absorption of light by a particle and it is 
smaller than 1 if some photons are absorbed by a 
particle. Clearly, the probability of photon survival, P, 
is equal to the single-scattering albedo (Chandrase- 
khar, 1950): 


(3.25) 


Osca 
wo = a (3.26) 
It should be pointed out that P= Gavbs/Gext is a small 
parameter for water clouds (see Fig. 5) in the visible 
and P=wo=1—f=1 in this case. The value 
P(Q)=x()/P is called the phase function. It is equal 
to one if the probability of light scattering does not 


depend on the angle (see Eq. (3.25)). Note that the 
value of x(Q) is sometimes also called the phase 
function or the scattering indicatrix (Ishimaru, 1978; 
Van de Hulst, 1980). 

Light scattering by water droplets is azimuthally 
symmetric due to their spherical shape, thus 


F p(0)sin0d0 = 1. (3.27) 
0 

The phase function p(0) (or the scattering indicatrix) 
can be calculated from Mie theory (Shifrin, 1951). It 
does not depend on the concentration of particles, by 
the definition, but on their refractive index and size. 
Functions p(0), obtained from Mie theory, for differ- 
ent particle-size distributions are presented in Fig. 7 at 
2=0.65 um. Main features of phase functions of water 
clouds in visible are (see Fig. 7): 


e sharp forward—backward asymmetry; 

e weak dependence on the size of droplets in the 
range of scattering angles 20—60°; 

e enhanced scattering near rainbow 0, (approximately 
138° for water droplets) and glory 0, (180°) scat- 
tering angles (Tricker, 1970; Greenler, 1980; Kon- 
nen, 1985; Spinhirne and Nakajima, 1994); 

e strong forward peak. 
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Fig. 7. Phase functions of water clouds obtained from Mie theory for the same conditions as in Fig. 3. Data for asr equal to 6 and 16 um are also 


given for the comparison. 


Amplitudes of peaks at angles 0=0, 0,, 2 can be 
used as indicators of the size of droplets (Van de 
Hulst, 1957; Shifrin and Tonna, 1993; Kokhanovsky 
and Zege, 1997b). 

Approximate equations for the phase functions of 
water clouds, derived using the geometrical optics 
concept, have been presented by many authors (e.g., 
see Shifrin, 1951; Van de Hulst, 1957; Nussenzveig, 
1992; Zege et al., 1993; Kokhanovsky and Zege, 
1997b; Grandy, 2000). A detailed comparison of the 
geometrical optics and Mie theory was performed by 
Liou and Hansen (1971), who conclude that the ray 
optics results have a high accuracy for large droplets 
with size parameters larger than 400. For some 
angular regions, errors are small even for compara- 
tively small droplets (see, e.g., the region 20—60° in 
Fig. 7). Geometrical optics equations are not very 
simple even for spherical droplets (Kokhanovsky, 
2001a). 

Thus, we found, parameterizing Mie theory results, 
that the phase function of a cloud can be presented 
approximately by the following equation: 


p(0) = Oexp(—C0) + X` biexp[—B,(0 — 6,)"), 
i=l 


where Q=17.7, C=3.9, 0 is given in radians and 
constants b;, P; 0; are given in Table 3. The short- 


coming of this equation is that it does not account for 
the influence of the size of droplets on the phase 
function (see Fig. 7). However, it can be useful for 
general modelling of single and multiple light scatter- 
ing in cloudy media (Zege et al., 1993, 1995; Katsev 
et al., 1998). 

The solution of the radiative transfer equation, 
which describes the radiative transfer through cloudy 
media, is simplified if one uses the expansion of the 
phase function in series: 


o0 


p(0) = S aP, (cos0). 


s=0 


(3.28) 


Here P,(cos 0) is the Legendre polynomial. However, 
note that in practice, the number of coefficients, a,, 
could be taken to be finite. It is approximately equal to 
2kd, where d is the average diameter of particles and k 
is the wave number (Kokhanovsky, 1997). The coef- 


Table 3 

Parameters b;, P; and 0; 

i b; Pi 0; 
1 1744.0 1200.0 0.0 
2 0.17 75.0 25 
3 0.30 4826.0 T 
4 0.20 50.0 T 
5 0.15 1:0 T 
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ficients, as, for cloudy media calculated with Mie 
theory at 4=0.65 um for the PSD (Eq. (2.3)) with 
u=6 are presented in Fig. 8. Approximate equation 
for as, based on the combination of the geometrical 
optics and the Fraunhofer diffraction, were obtained 
by Kokhanovsky (1997, 2001a). They can be used to 
avoid long numerical calculations with Mie theory as 
Xef— 0. 

We emphasize that both geometrical optics results 
(Kokhanovsky, 1997) and data presented in Fig. 8 
show that the number s"™*, which corresponds to 
the largest coefficient aj"*“, is equal approximately 
to 4xXe9/5. Also we have approximately s™*a}"*“=2. 

It follows from Eq. (3.28) using the orthogonality 
of Legandre polynomials 


_2s+1 
ay = 5 


| p(0)Ps(cos0)sinOdd. (3.29) 
0 
In particular, we have aọ=1 (see Eq. (3.27)) and 
a, =3g, where 


1 [* : 
g= 5S p(0)sin@cosddé (3.30) 
0 
is the so-called asymmetry parameter, giving the 
average cosine of the photon-scattering angle. Note 
that the value of g determines (together with the 


extinction coefficient) the coefficient of photon dif- 
fusion, D=[3¢.x:(1 — g)]~ l ina cloudy medium. 

We stress that it is the value of integral (Eq. 
(3.30)) (and not the phase function itself), which 
determines a number of important radiative charac- 
teristics of clouds (Zege et al., 1991). For instance, 
we have for the value of the total cloud reflectance 
or the spherical albedo in the case of a semiinfinite 
cloud (Zege et al., 1991; Kokhanovsky, 2001a): 
ro = exp(—4,/$/3(1 — g)). Thus, the asymmetry 
parameter needs a special attention. We will consider 
this parameter in the next section in more detail. 


3.1.4. Asymmetry parameter 

Let us represent the phase function p(0) of a large 
particle (xef — œ, 2Xef(n— 1)—> œ) as (Kokhanov- 
sky, 2001a) 


a? P” (0) + o&, p° (0) 
oD 4 g6 : 


sca T “sca 


p(0) = (3.31) 


Here p”(0) is the phase function associated with the 
Fraunhofer diffraction of light by a droplet; p%(0) is 
the phase function associated with rays, which pene- 
trate the particle, reflect, and refract inside the droplet; 
o&, and a2, are scattering coefficients associated with 
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S 


Fig. 8. Coefficients a, obtained for the same conditions as in Fig. 3. Data for aep equal to 6 um are also shown. 
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these processes. Note that Eq. (3.31) corresponds to 
the Van de Hulst’s (1957) localization principle. It 
holds only approximately and ignores the possible 
interference of waves (Glautshing and Chen, 1981) 
that originate due to separate diffraction and geo- 
metrical optics scattering processes. 

It follows from Eqs. (3.30) and (3.31) 


D -D G „G 
FscaS T Osca 


e= eron = 
where 
1 T 
P=5 f p?(0)sindcos0dd, (3.33) 
0 
=i f p%(0)sindcos0dd, (3.34) 
0 


One can obtain (Van de Hulst, 1957; Kokhanovsky, 
2001a,b) that g? = 1. Note that the function p%(6) 
does not depend on the size of large nonabsorbing 
particles (Van de Hulst, 1957; Kokhanovsky and 
Nakajima, 1998). Thus, g depends only on the 
refractive index, n, of water in this case. The value 
of n does not change considerably in the visible region 


of the electromagnetic spectrum (see Fig. 1). We will 
assume that n=1.333. Then it follows that g®= 
0.7686 (Kokhanovsky, 2001a) and g=0.8843 accord- 
ing to Eq. (3.32), where we accounted for the equality 
o2,=08%,, which also holds in the visible region of 
the electromagnetic spectrum, where absorption can 
be neglected (Van de Hulst, 1957). The value of 
asymmetry parameter g of water clouds, in practice, 
weakly depends on the size of particles and wave- 
length. This dependence can be approximated by the 
following equation (Kokhanovsky, 2001a): 


C 
Xef 


where C is the constant, which does not depend 
on the size of droplets (but possibly depends on 
their shape) and go = g(%er— œ). The integral for 
go is given by Van de Hulst (1957). It depends 
only on the refractive index of a particle and equal 
approximately to 0.8843 at n=1.333, as stated 
above. 

Results of calculations of the value C=(e0- 2x} a 
with Mie theory are presented in Fig. 9 for the 
particle-size distribution (Eq. (2.3)) at u=6, aef=4 
um, and 2 <1 um. It follows that C = 0.5. Thus, we 
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Fig. 9. Spectral dependence C(A) obtained for the same conditions as in Fig. 3. 
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have from Eq. (3.35) the value of G=1— g, which is 
often called the co-asymmetry parameter: 


1 
G = 0.12 +—,. 
22h 


(3.36) 


It is the value of G and not g itself which plays an 
important role in the radiative transport in cloudy 
media. For instance, values of g equal to 0.75 (typical 
for ice clouds) and 0.85 (typical for water clouds) are 
comparatively close to each other, but corresponding 
values of G (0.25 and 0.15, respectively) differ. This is 
a noteworthy feature for comparing cloudy media 
with different values of asymmetry parameters. 

It is important to understand the physical meaning 
of the parameter G. For this, we use the expansion 
cos0= 1 — (67/2)+... in Eq. (3.30). Then it follows 
that G= (0°)/2, where (0°) =} fọ 0°p(0)sinOd0 is the 
averaged square of the scattering angle. We neglect 
high-order terms in the expansion of cos@ which is 
possible due to a peaked character of the phase 
function in Eq. (3.30) (see Fig. 7). 

The average value of the squared scattering angle, 
therefore, governs the radiative transport in cloudy 
media. In particular, semiinfinite clouds with close 
values of ratios f/(0°) also have close values of the 
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total reflectance. Smaller values of (0°) mean that 
photons penetrate to larger depths in cloudy media 
before their escape back to the free space, where the 
source of radiation is located. This also makes the 
total path lengths of photons in cloudy media larger, 
which in turn produces the larger total absorption of 
radiation inside clouds with larger values of (07). 

The total fraction of radiation absorbed by a semi- 
infinite cloud is given by the following expression 
(Kokhanovsky, 2001a): 


A=1 oof 4 E 
or 
ai 
3(0°) 
as f — 0. Clearly, the ratio 
_ A 
i =y 


can be used as an estimation of the average number of 
photon-scattering events in a cloud layer. In particular, 
' 2 
@2 


we get n — 4, /5 as B — 0. Smaller values of p 
and (0°) give larger values of 7 as one might expect. 
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Fig. 10. Asymmetry parameter obtained for the same conditions as in 
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Fig. 3. Results for the effective radius a.¢=6 um are also shown. 
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Fig. 11. Error of the geometrical optics approximation for the asymmetry parameter obtained from data presented in Fig. 10 for effective radii 4 


and 6 um. 


Eq. (3.36) can be used only for nonabsorbing 
channels (e.g., in visible). However, it can be modi- 
fied in a full analogy with Eq. (3.23) to account for 
light absorption at 2> 1 um (Kokhanovsky, 2001a): 


1 
G= 0.12 + —s5 — taer. (3.37) 
2. 


Xef 


We obtained fitting Eq. (3.37) to Mie theory results 
that e=0.15. Comparisons of calculations with simple 
Eq. (3.37) and Mie theory are presented in Fig. 10. 
They were performed for the particle-size distribution 
(Eq. (2.3)) at w=6, dep=4 um, and age=6 um for 
wavelengths, 2 < 2.4 um. It follows that the accuracy 
of Eq. (3.37) (see Fig. 11) is better than 5% at a.p=6 
um in visible and near infrared for å < 2.3 um. It is 
also better than 5% for droplets with aep=4 um for 
values of A < 2.0 um. 

Another important characteristic of the phase func- 
tion of a cloudy medium is the probability of light 
scattering in the backward hemisphere: 


1 T 
F= J v@sindcosoao, 


It was shown (Kokhanovsky et al., 1998a) that it holds 
approximately for water clouds: 


1 


F = 0.03 +—,. 
5x23 


It suggests that approximately 97% of light is 
scattered by a local volume of a cloudy medium in 
the range of angles smaller than 90°. The backscatter- 
ing signal is, therefore, quite low. It increases, how- 
ever, due to multiple-scattering processes that take 
place inside a cloud (Bissonnete et al., 1995; Zege et 
al., 1995, 1998a,b). 


3.2. Ice clouds 


3.2.1. Extinction coefficient 

Local optical characteristics of ice and mixed 
clouds cannot be calculated so easily as for water 
clouds. In particular, the Mie theory (Shifrin, 1951) 
cannot be used due to the complex shape and internal 
structure of ice crystals (Takano and Liou, 1989; 
Macke, 1993, 1994; Macke et al., 1996, 1998; Mitch- 
ell and Arnott, 1994; Yang and Liou, 1998; Yang et 
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al., 2000, 2001). Main results obtained in the optics of 
crystalline clouds were summarized by Volkovitsky et 
al. (1984) and Liou (1992, 2002). 

Ice crystals are usually much larger than the wave- 
length of the incident radiation. Thus, the extinction 
cross section, C.x;, does not depend on the wavelength 
and the refractive index of particles (Shifrin, 1951; 
Van de Hulst, 1957): 


Cext = 25, (3.38) 
where s is the cross section of a particle, projected on 
the plane perpendicular to the incidence direction. It 
follows for N identical crystals in a fixed orientation 
that 


ext = 2Ns, (3.39) 
where N is number concentration of particles and ox; is 
the extinction coefficient. However, crystals differ by 
their shape, size, and orientation; thus, the extinction 
coefficient can be calculated as 


Gext = 2N5, (3.40) 


where 5 is the average cross section of particles. This 
equation can be written as 


Oext = 24>; (3.41) 


NS 


where Cy is the volumetric concentration of crystals 
and V is the average volume of crystals (C, = ND). It 
follows for convex crystals of the same shape in 
random orientation (Van de Hulst, 1957) that 5= 4s, 
where X is the average surface area of particles. There- 
fore, we get: 


3Cy 
at=; 42 
Tezi 2aef (3 ) 
where 

3V 
T=- 3.43 
det = (3.43) 


is the effective radius of particles. This result is similar 
to Eq. (3.7). Only a portion of ice crystals (e.g., 
hexagonal plates and columns) are convex. For this 
and other reasons, one should expect that Eq. (3.42) 


will provide us with only a coarse approximation. It 
follows from Eq. (3.42) that 


(3.44) 


Eq. (3.44) suggests a method for indirect determi- 
nation of the effective radius of particles in complex 
crystalline cloudy media. 


3.2.2. Absorption coefficient 

The absorption cross section of a single crystal at 
small absorption (&Xaef— 0) can be found in a full 
anology with Eq. (3.16): 


Cabs = Bal, (3.45) 
where V is the volume of a crystal, «=4zy// and B is 
the coefficient proportionality, which depends on the 
shape and the real part of the refractive index of 
particles (Kokhanovsky and Macke, 1997), but not 
on their size (see Eq. (3.14)). 

For an ensemble of crystals of identical shapes: 


Cabs = BaC, (3.46) 


and for N distinct shapes of crystals in a cloud: 


N 
Tabs = 2 > BCy. (3.47) 
j=l 

Values of B; do not depend on the imaginary part of 
the refractive index and the size of crystals, but only 
on the real part of the refractive index, n, of crystals 
and their shapes. The most probable values of B= 
(1/C,)5*_ ;B,C,, in crystalline clouds can be found 
from experimental measurements of ratios o4p,/Cyo in 
ice clouds of different microstructures. 

The probability of photon absorption f= dap5/Cext 
can be found from Eqs. (3.42) and (3.47): 


B = Sader, (3.48) 
where 
2 A 
a B;Cyj (3.49) 
3G 2 


Using the model of fractal particles (Macke, 1994), 
one can obtain that £ = 1.8. Due to the complexity of 
crystal shapes in ice clouds, the accurate value of 2 
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cannot be found theoretically. Thus, measurements of 
this parameter are needed. 

It is readily apparent from Eq. (3.48) that the 
single-scattering albedo 
wo = 1 — Haag (3.50) 
for crystalline media. Note that the formulae presented 
here are valid only as a&age— 0. Generalization of 
these equations for arbitrary light absorption by drop- 
lets was given by Kokhanovsky and Macke (1997). 
Most recent parameterizations were developed by 
Yang et al. (2000) and McFarquhar et al. (2002). 


3.2.3. Phase function 

The phase function of ice clouds in the visible is 
the average over ensemble of phase functions of 
crystals with different shapes (Takano and Liou, 
1995). It depends on the size of crystals in the 
small-angle scattering region and primarily on shape 
and structure of crystals at larger scattering angles. 
The phase function of ice clouds cannot be modelled 
(Liou and Takano, 1994) by spherical polydisper- 
sions. This is in contrast with integral light-scattering 
and absorption characteristics of crystalline clouds, 


phase function 
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which in many cases (see, e.g., Eqs. (3.21) and (3.48)) 
do can be modelled by a collection of spherical 
particles (Grenfell and Warren, 1999). 

The phase function of hexagonal cylindrical crys- 
tals outside the diffraction region, presented in Fig. 
12, was obtained using the Monte-Carlo ray tracing 
approach (Macke, 1994) for identical hexagonal ice 
cylinders of length 0.5 mm in random orientation at 
wavelength 0.5 um. The side of the cross section was 
taken to be equal to 0.08 mm. Dominant features of 
the phase functions of this type are enhanced scatter- 
ing at angles of 22° and 46°, corresponding to the so- 
called halos. Greenler (1980) states that he observes 
22° halos in Wisconsin on 70 to 80 days of a typical 
year. The second halo at 46° is of a rare occurrence in 
natural conditions due to the presence of crystals of 
other shapes, which obscure halo phenomena. Both 
rainbow and glory (Tricker, 1970; Greenler, 1980; 
Konnen, 1985) are absent for ice clouds. 

The backward cross section of hexagonal ice 
crystals of arbitrary orientation is calculated for visi- 
ble light by means of a ray-tracing code by Borovoi et 
al. (2000). It is shown that backscattering of the tilted 
crystals is caused by a corner reflectorlike effect. A 
very large peak of backscattering is found for a tilt 


hexagonal cylinder 
fractal particle 


100 120 140 160 180 
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Fig. 12. Phase functions of hexagonal ice cylinders with the aspect ratio(length/size of the side of the hexagonal cross section) of 5.88 and fractal 
particles in random orientation at the wavelength 0.5 um. Only the geometrical optics contribution of both phase functions is shown. 


210 A. Kokhanovsky / Earth-Science Reviews 64 (2004) 189-241 


angle of 32.5° between the principal particle axis and 
the incidence direction. This peak is caused by multi- 
ple total internal reflections for part of the rays that are 
incident upon the skewed rectangular faces. 

Realistic phase functions of a crystalline cloud that 
account also for diffraction of light at small angles and 
size/shape distribution of crystals are tabulated by 
Takano and Liou (1989) and Liou (1992) and are 
characterized by an asymmetry parameter equal to 
0.75. 

Calculation of phase functions of crystals of differ- 
ent shapes can be avoided by introducing the single 
fictive particle with the phase function similar to the 
phase function of an ensemble (Kokhanovsky, 2001a, 
2003a,b,c). Particles with random stochastic surfaces 
were studied by Peltoniemi (1993) and Muinonen et 
al. (1996). The model of a fractal particle was intro- 
duced by Macke and Tzschichholz (1992), Macke 
(1994), and Macke et al. (1996). Depending on the 
parameters of roughness, these models can quite well 
describe the phase functions of ice crystals (Macke, 
1994; Kokhanovsky, 2003a,b,c) and the mineral frac- 
tion of the atmospheric aerosol (Volten et al., 2001; 
Kokhanovsky, 2003a,b,c). 

Phase function of a fractal particle, which may be a 
better description of complex shapes in crystalline 
clouds, is presented (along with that of a randomly 
oriented hexagonal cylinder) in Fig. 12. Data were 
obtained using the Monte-Carlo ray-tracing code 
(Macke, 1994; Macke et al., 1996). The fractal phase 
function is characterized almost the same value of the 
asymmetry parameter (g = 0.74) as the phase func- 
tion of randomly oriented hexagonal crystals given in 
the same figure. 

Linear combinations of functions presented in Figs. 
7 and 12 with weights depending on the concentration 
of spherical particles, hexagonal cylinders and irreg- 
ularly shaped particles in the cloud can be used to 
describe natural clouds. Note that other regular 
shapes, such as plates, can also contribute to the total 
phase function of an elementary volume of a cloudy 
medium. Measurements of phase functions of crystal- 
line clouds were presented by Volkovitsky et al. 
(1984), Gayet (1998), Saunders et al. (1998), Barkey 
et al. (1999), Baran et al. (2001) and Jourdan et al. 
(2003). 

Horizontal orientation of crystals may complicate 
the phase function in up to 40% of clouds (Chepfer et 


al. (1999)). In particular, the phase function of hori- 
zontally oriented crystals (e.g., plates) depends on 
directions of light incidence and observation sepa- 
rately and not on a single parameter (scattering angle) 
as for spherical or randomly oriented particles. Such 
clouds belong to a broad class of the so-called 
anisotropic media. The radiative transfer in ice clouds 
with horizontally oriented crystals was studied by 
Liou (2002). This is a major topic of modern cloud 
optics research. Note that clouds with horizontally 
oriented crystals are characterized by dichroism and 
birefringence (Van de Hulst, 1957). 


3.2.4. Asymmetry parameter 

Asymmetry parameters of crystals depend on the 
shape, but not on their size in the visible (Macke et al., 
1996; Kokhanovsky and Nakajima, 1998), due to 
large size of ice crystals. Size-dependent wave cor- 
rections (see, e.g., Baran and Havemann, 1999), such 
as given by the second term in Eq. (3.36), can be 
neglected for crystals as 4 — 0. 

The real part of the refractive index of ice crystals 
is also important, but it varies only slightly in the 
visible (see Fig. 1); this dependence can be neglected. 
The asymmetry parameter for ice clouds will be the 
average value over an ensemble of particles of differ- 
ent shapes. Macke et al. (1998) found from extensive 
numerical calculations that the asymmetry parameter, 
g, is in the range 0.79—0.85 for columns, 0.83—0.94 
for plates, and 0.74 for polycrystals, represented by a 
fictive fractal particle in the visible. It is in the range 
0.83—0.87 for spherical droplets. We see, therefore, 
that the value of g for cloudy media with ice crystals 
should be somewhere between 0.74 and 0.87. Larger 
values of g indicate that ice clouds are composed of 
plates only, which is never the case in natural con- 
ditions. We see that the co-asymmetry parameter, 
G=1-—g, changes from 0.13 for water clouds to 
0.26 for irregularly shaped particles. Actual values 
of g for a given cloud can be found only from direct 
measurements in natural ice clouds such as were 
performed by Garett et al. (2001), who found that 
the value of g is in the range 0.73—0.76, depending on 
the cloud area under investigation. They obtained 
values of g from measurements of the phase function 
inside the ice cloud in the angle range 10—175°. In 
one case, the value of g reached 0.81, but sublimation 
was noted, rendering it to be a typical case. The 
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asymmetry parameter of ice clouds does not vary 
considerably for visible light and, on average, is equal 
to 0.745, a value close to that for a fractal fictive 
particle and to the value for hexagonal crystals, which 
give the phase function, presented in Fig. 12. 

We suggest that an asymmetry parameter value of 
0.74 should be used for theoretical modelling of light 
propagation in crystal clouds. The corresponding 
value of G is equal to 0.26, which does not depend 
on the size of crystals because they are large as 
compared to the wavelength in the visible range of 
the electromagnetic spectrum. G does not depend on 
the shape of crystals because for a given cloud, we 
have a statistical and very broad distribution of 
shapes, which produce in the end a saturated value 
of the asymmetry parameter for a completely chaotic 
scattering. This scattering can be modelled by a single 
fractal fictive particle (Macke, 1994; Kokhanovsky, 
2003a). It should be stressed, however, that the 
calculation of g involves also averaging over scatter- 
ing angle (see Eq. (3.30)). Thus, the model of chaotic 
scattering is more appropriate for the value of g than 
for the phase function itself. The phase function 
should be modelled as a combination of light scatter- 
ing by regular and irregularly shaped particles with 
different weights, as discussed above (see, e.g., Liou, 
2002). Again, such a model can be established only 
from measurements performed in natural clouds. 

Now we consider briefly the case of mixed ice— 
water clouds for which the value of G is 


Gnt (3.51) 


where G; and Gw are co-asymmetry parameters of ice 
and water particles, respectively, and T is the fraction 
of water droplets in a mixed cloud, defined as the ratio 
of number of droplets to a total number of particles in 
a mixed cloud. G; and Gw do not vary significantly in 
mixed clouds; thus, in the first approximation, we can 


assume constant values for them. Measurements by 
Garett et al. (2001) suggest that G;=0.26 and 
Gy =0.13. Then we have from Eq. (3.51): 


G 


jo 
0.13 


2 (3.52) 


or approximately Y=(4(07))—2. The correlation 
between values of G and Y was confirmed experimen- 
tally (Garett et al., 2001). Therefore, measurements of G 
(or, alternatively, (0°)) in mixed clouds can yield the 
value of a liquid fraction, Y. 

It is evident (Kokhanovsky and Macke, 1997; 
Kokhanovsky, 2001a) that in the near-infrared range 
of the electromagnetic spectrum the larger size of 
crystals generally leads to larger asymmetry parame- 
ter. Detailed calculations of optical characteristics of 
ice crystals in the infrared for different shapes of 
crystals were performed by Zhang and Xu (1995) 
and Macke et al. (1998). Parameterization schemes for 
local optical characteristics of ice clouds were devel- 
oped by Ebert and Curry (1992), Fu and Liou (1993), 
Fu (1996), Kokhanovsky and Macke (1997), Yang et 
al. (2000), Harrington and Olsson (2001), and McFar- 
quhar et al. (2002). Selected parameterizations of 
cloud local optical characteristics given in this review 
are summarized in Table 4. 


4. Radiative characteristics 
4.1. The visible range 


Let us consider now the global optical character- 
istics of clouds such as their reflection and trans- 
mission functions (Van de Hulst, 1980). They can be 
measured remotely by airborne, satellite and ground- 
based radiometers and spectrometers (Kondratyev and 
Binenko, 1984; Liou, 2002). The aim of this section is 


Table 4 

Integral local optical characteristics of clouds in visible and near infrared 

LOC Water cloud Ice cloud Range of validity 

Oabs (5/3)a(1 — odep)Cy Bac, dep — 0 

ext (3/2depC(1 + 1.x 27) (3/2dex)Cy Xep— 00, 2Xer(n — 1)— o 

l-g 0.12 +0.5xe¢ 77 — 0.150der 0.26 Xet — 00, 2Xer(n — 1) 00, «aet — 0 


Here #=42y/A and xef = 2na A. The parameter B depends on the shape of crystals in a cloud (Kokhanovsky and Macke, 1997). It follows for 
fractals: Bx 1.8. 
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to introduce simple analytical formulae, which can be 
used to relate global optical characteristics of clouds 
to their local optical characteristics discussed in the 
previous section. This will also mean that we obtain 
direct relationships between observable quantities, 
such as cloud reflection and transmission functions, 
with microphysical parameters of clouds (e.g., the size 
and refractive index of particles, see Table 4). 

We start from the study of the reflection function of 
a cloudy medium under an assumption that the cloud 
can be represented as a homogeneous plane-parallel 
layer. The absorption of light by particles is neglected 
in this section. This is a standard assumption in the 
visible range of the electromagnetic spectrum. How- 
ever, in some cases the absorption of light by water 
clouds in visible cannot be neglected. This is due to 
various sources of pollution (Melnikova and Mikhai- 
lov, 1994, 2000; Asano et al., 2002; Vasilyev and 
Melnikova, 2002). 

The reflection function of a cloud R(Vo,V,@) is 
defined as the ratio of reflected light intensity 
L(V, p) for the case of a cloud to that of an ideal 
Lambertian ideally white reflector 


= I(00, 0, o) 


R(Io, 0 ee . 

( 0; 9) I* (0o) ’ (4 1) 
where 
I* (0o) = Fosvo (4.2) 


is the intensity of light reflected from the ideally white 
Lambertian reflector, zF is the solar flux on the area 
perpendicular to the direction of incidence, Wo is the 
solar angle, 7 is the observation angle, and ¢ is the 
relative azimuth between solar and observation direc- 
tions (Minin, 1988). Also, we will use u= | cosv |, 
Lo =cosvo. For the Lambertian ideally white reflector 
from Eq. (4.1), R = 1, which, by definition, does not 
depend on the viewing and observation geometry. 
Although clouds are white when looking from 
space, their reflection function R(v,0,@) is not equal 
to one, but depends on the viewing and observation 
geometry. Results of calculations of the reflection 
function of an idealized semiinfinite nonabsorbing 
water cloud R°, (Vo, 0, p) at the wavelength, 1=0.65 
um and the nadir observation are presented in Fig. 13. 
Calculations were performed for the gamma particle- 
size distribution, given by Eq. (2.3) at u=6. Values of 


effective radius in Fig. 13 were 6 and 16 um, which 
cover the typical range for natural water clouds. 
R(Vo,0,@) can be smaller and larger than 1, depending 
on the incidence angle, which implies that depending 
on particular viewing geometries, clouds can be even 
more reflective than the ideally white Lambertian 
surface. This is mostly due to peculiarities of the 
phase function of cloudy media [e.g., in the back- 
scattering (V = vop = m) region]. It follows from 
Fig. 13 that in the range of solar angles 30—60° and 
nadir observation, the reflection function of a semi- 
infinite water cloud is almost equal to the reflection 
function of an ideally white Lambertian reflector. It 
differs from 1 not more than by 10% for these geo- 
metries (Zege et al., 1991). 

The reflection function R2(o,0,) can be repre- 
sented by the following simple approximate equation 
(Kokhanovsky, 2002a): 


bı + bzcosðcosvo + p(0) 
~ (Vo V, p) 4(cosð + cosvo) 


(4.5) 


where 0 = arccos(— cos¥cosvo + sinvsinvocoso) is the 
scattering angle, p(0) is the phase function of a cloudy 
medium, and b; and bp are constants. Eq. (4.5) obeys 
to the reciprocity principle [R°.(9,0,2) =R°(0,00,0), 
see, e.g., Zege et al., 1991]. 

We found that b; = 1.48, b2 =7.76 for nadir obser- 
vations (Kokhanovsky, 2002a). Comparison of the 
approximation and exact data in Fig. 13 shows that 
the accuracy of Eq. (4.5) is better than 2% at Vo <85°, 
0=0°. Constants b; and bz for other viewing geo- 
metries can be found using parametrizations of results 
obtained from the exact radiative transfer codes (see, 
e.g., Mishchenko et al., 1999). Approximations for the 
function R°(9,0,¢~) at V= 0° were also obtained by 
Melnikova et al. (2000a,b) and Kokhanovsky (2003c). 

Eq. (4.5) can be used to find the reflection function 
of a finite cloud R(u,0,~,T) via the following formula 
(Germogenova, 1963; Van de Hulst, 1980, Minin, 
1988): 


R(u, Hos P, 1) = R? (u, uo, P) — t(t)Ko(u)Ko( uo), 


4.6 
where (4.6) 


1 
0.75t(1—g) +4 
is the global transmittance of a cloud, g is the asym- 
metry parameter studied above in detail, and Ko() is 


t (4.7) 
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Fig. 13. Reflection function of an idealized semiinfinite nonabsorbing cloud R°.(0,Jo,0) obtained from the exact radiative transfer code 
(Mishchenko et al., 1999) and approximation (4.5) at the wavelength 4=0.65 um and the effective radii of droplets a.¢=6 and 16 um. It is 
assumed that particles in a cloud are characterized by the gamma particle-size distribution (Eq. (2.3)) with the parameter = 6. 


the escape function (Van de Hulst, 1980). Eqs. (4.6) 
and (4.7) are valid only in the absence of light 
absorption in a cloud. The escape function is defined 
via the solution of the characteristic integral equation 
(Van de Hulst, 1980). Note that the parameter « in Eq. 
(4.7) is defined as follows (Sobolev, 1972): 


a= 3 | Ko(u)pe du. (4.8) 


The escape function, Ko(), only weakly depends on 
the cloud microstructure and can be presented by the 
following equation (Zege et al., 1991; Kokhanovsky, 
2001a): 


Ko(u) ==(1 + 2u). (4.9) 


3 
7 
The function Ko(t9) calculated with exact radiative 


transfer code with g equal to 0.75, 0.85, and 0.9 in the 
case of Heney—Greenstein phase function 


p(0) = pal 1)g'P;(cos0), (4.10) 


where P} is the Legendre polynomial and, according to 
Eq. (4.10), is presented in Fig. 14. We see that Ko() 
almost does not depend on g at u > 0.2 (0 <78°). This 
is the case even for g=0 and for the Mie-type phase 
functions (Kokhanovsky, 2001la). There is some 
dependence of the escape function on the microstruc- 
ture of the cloud at the range of observation angles 
V =80-—90°, but the cloud-top inhomogeneity plays a 
role at such grazing observation angles, preventing use 
of the idealized plane-parallel layer approximation. 

The variability of Ko(u) at ~=0.2—1.0 for different 
values of g=0.75—0.9 is well inside the 2% sensitiv- 
ity range. This coincides with the error of Eq. (4.9) at 
= 0.2, being smaller than 2%. 

Note that function (4.9) also describes the angular 
distribution of solar light transmitted by a cloud (Zege 
et al., 1991; Kokhanovsky, 2001a). 

The substitution of Eq. (4.9) into Eq. (4.8) yields: 


15 
= — ~1.07 4.11 
“=T (a 
independently of cloud microstructure. The value of «, 
numerically calculated by King (1987), assuming the 
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escape function 
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Fig. 14. Escape function, calculated with exact radiative transfer code for the Heyney—Greenstein phase function at g=0.75, 0.8 and 0.9 


(Yanovitskij, 1997) and with approximation given by Eq. (4.9). 


fair weather cumulus cloud model, is given approx- 
imately by 1.07 in agreement with our estimation. 
King (1987) used Mie theory to find the phase 
function of a cloudy medium for the calculation of 
Ko(o) and « via Eq. (4.8). Yanovitskij (1997) found 
the same value of « for Heney—Greenstein phase 
functions with asymmetry parameters in the range 
0.0—0.9. This supports the approximation of a fixed 
value of « in Eq. (4.7), given by Eq. (4.11), inde- 
pendently of cloud microphysical parameters. 

The error of Eq. (4.6) is smaller than 1% at t > 10 
(King, 1987). It rapidly increases with t. Kokhanov- 
sky and Rozanov (2003) proposed to modify Eq. (4.6) 
to diminish the error of this equation for smaller q. In 
particular, they proposed to use instead of ¢ in Eq. 
(4.6): 


= —7, 7 = (4.86 — 13.08 upo + 12.767 u$) /t. 


This modification allows to increase the accuracy 
of approximate formulae considered for smaller t 
(down to t=5). This range of optical thicknesses is 
that most frequently observed in water clouds both 
with satellite and ground-based techniques (see Fig. 


15, prepared from data given by Trishchenko et al., 
2001). Kokhanovsky and Rozanov (2003) found that 
the error of the modified reflection function at the 
nadir observation is smaller than 5% for solar angles 
smaller than 74°. 

Eq. (4.6) is readily modified to account for the 
Lambertian light reflection from the underlying sur- 
face (Sobolev, 1972): 


, rTG)T (m) 
Ru, to, 0,1) = R(H, po, Qt) + SG 


(4.12) 


where R is the reflection function of a Lambertian 
surface-cloud system, rs is the spherical albedo of the 
Lambertian surface, which may depend on the wave- 
length, R = R (r,=0), 


T (u) = Kolu) 


is the diffuse transmittance of a cloud layer (Sobolev, 
1972) and r=1—t is the spherical albedo of a non- 
absorbing cloud. 

We have neglected the direct solar light term (Liou, 
2002) in Eq. (4.12) due to a large thickness of clouds 
under consideration. 


(4.13) 
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Fig. 15. Frequency of registration of different optical thicknesses of cloudy media obtained from satellite and ground measurements as presented 


by Trishchenko et al. (2001). 


Finally, substituting Eqs. (4.6) and (4.13) into 
Eq. (4.12), it follows for the reflection function of 
a Lambertian surface-nonabsorbing cloud system 
(Kokhanovsky et al., 2003): 


R(u, Hos P, T) = RY (u, Mo, P) 


7 E Kal Kolt) 


(4.14) 


This formula can be used as a basis for the semi- 
analytical cloud retrieval algorithm (Kokhanovsky et 
al., 2003). Note that (u, Ho,0,1) = R? (u, 40,9) at 
r=1. 


4.2. The near-infrared range 


Unfortunately, the relatively simple Eq. (4.6) can- 
not be used to calculate the reflection function of a 
finite cloud in the near-infrared region of the electro- 
magnetic spectrum due to the presence of absorption 
bands of liquid water (see Fig. 2). Alternatively, the 
following formula applies (Germogenova, 1963; Van 


de Hulst, 1980; King, 1987; Nakajima and King, 
1992; Vasilyev and Melnikova, 2002): 


R(u, Ho, P, T) — Reali Ho, o) 


mle~?"" 


-q pear KHK), (4.15) 


where y is the so-called diffusion exponent, K(u) and 
R „~ are the escape function and the reflection function 
of an absorbing semiinfinite medium with the same 
local optical characteristics as a finite layer under 
study, respectively. Eq. (4.15) accounts for the influ- 
ence of light absorption on the reflection function of 
clouds. The reflection function decreases if additional 
absorbers are present in cloud droplets. Constants / 
and m are defined by the following integrals (Van de 
Hulst, 1980): 


1 
l= 2 | i (n)ndn, (4.16) 
-1 


m=2 f K(n)i(—n)ndn, (4.17) 
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i(n) being the angular distribution of light in deep 
layers of a cloud, defined as the so-called asymptotic 
regime (Sobolev, 1972). 

Functions R (4, Ho,), K(u), and constants m and / 
have the following asymptotic forms when light 
absorption by droplets is relatively small (wg — 1) 
(Van de Hulst, 1980): 


R° (u, Ho, P p) 


l-o 
-4r Kololo), (4.18) 


Ro (H, ko, P 


1— 
K() = nwi- nOA ee At 2 5) (4.19) 
T ea! (4.20) 
m= ; 
3(1-g) 
E T —™ (4.21) 
3(1— g)’ l 
and y 3(1 — wo)(1 — g) as the single-scattering 


albedo, wp—1. Thus, at œo=1, R..=R°,, K(w) = Ko(o, 
m=1=0, and Eq. (4.15) transforms into Eq. (4.6). 
Note that Eqs. (4.18)—(4.21) follow from the 
asymptotic analysis of the radiative transfer equation 
(Van de Hulst, 1980). Integration of Eq. (4.18) with 
respect to all angles yields r,,=1—y, where y= 4 


1—ao 
EEA) and ry% 


is the spherical albedo of a weakly 


absorbing semiinfinite cloud (Van de Hulst, 1980). 
Thus, the parameter y can be interpreted as a fraction 
of photons absorbed in a weakly absorbing semi- 
infinite cloud. It depends both on the single-scattering 
albedo and the asymmetry parameter. Larger values of 
g imply that photon scattering increases at small 
angles. Thus, the photon path length before its escape 
from the medium is also increased. As a consequence, 
this results in increased light absorption in the 
medium. Clouds having larger values of g, therefore, 
absorb more light. Media having different values of 
œ and g, but the same values of y, have the same 


1—00 


og is called 


values of r. The parameter s = 


the similarity parameter in the radiative transfer text- 
books (Van de Hulst, 1980). It is a useful parameter, 
describing the optical properties of clouds. We see that 
y—4s as @wo—> 1. Substituting y into Eqs. (4.18)— 
(4.21) yields: 


Roo (H, Ho, P) = R? (H, Ho, P) — YKo(u)Ko (uo), 
(4.22) 

Ku) = Ko (1- 2), (4.23) 

m = 2y, (4.24) 

and 

l= 1- o. (4.25) 


Eqs. (4.22)-(4.25) were derived assuming that 
œo —> 1. Alternatively, the right-hand sides of Eqs. 
(4.22)—(4.25) give us the first terms of the expansion 
of corresponding functions with respect to y (Minin, 
1988). The accuracy of these equations decreases with 
y quite rapidly. The next terms of the expansions 
(~ y’) are given by Minin (1988) and Melnikova et 
al. (2000b). Still higher terms are not known or quite 
complex (Minin, 1988), but it has been shown that 
the following equations account approximately for 
higher-order terms (Rozenberg, 1962; Zege et al., 
1991; Kokhanovsky et al., 1998b, 2003): 


Ra (H, Ho, P) = Ro,exp[—vu(u, Mo, 9)], (4.26) 
mK (1) K (to) = [1 — exp(—2y)]Ko(“)Ko(Ho); 

4.27) 
l = exp(—ay), (4.28) 
where a viewing function u is defined by 
u(H, Ho, P) = eo (4.29) 


and does not depend on wp and t. The viewing 
function has only a weak dependence on the micro- 
physical characteristics of clouds (e.g., the droplet-size 
distribution). This follows from the low sensitivity of 
the functions R? (u,40,9) and Kolu) in Eq. (4.29) to the 
microphysical characteristics of clouds (see Figs. 13 
and 14). The physical sense of the function (4.29) has 
been discussed by Kokhanovsky (2002b). The accu- 
racy of Eq. (4.26) was studied by Kokhanovsky 
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(2002c). It is better than 10% at y < 0.8 for the nadir 
observation and solar angles in the range 0—75°. 

Eqs. (4.26)—(4.28) avoid the solution of integral 
equations (Van de Hulst, 1980; Yanovitskij, 1997) for 
the determination of functions R „(4,Ho,0) and K(x) in 
Eq. (4.15). They transform into exact asymptotic 
results (Eqs. (4.22)—(4.25)) as y— 0. 

Substituting Eqs. (4.27) and (4.28) into Eq. (4.15), 
we have: 


R(u, Hos P, T) = Ro (u, Ho, o) ig te ™Ko(u)Ko(uo), 
(4.30) 


where x=yt and the global transmittance, ¢, is given 
by (Kokhanovsky et al., 2003): 

sinhy 
on sinh(«y + x) net) 
Eq. (4.31) yields Eq. (4.7) at @p=1. The value of 
exp(— x) describes light attenuation in deep layers of a 
cloud (Van de Hulst, 1980). 

Eq. (4.30) was first proposed by Rozenberg (1962, 
1967), but he assumed that «=1 in Eq. (4.31), which 
is not consistent with the exact asymptotic result given 
by Eq. (4.7). 

The range of applicability of Eq. (4.26) with 
respect to higher values of y can be readily extended 
using the following correction (Kokhanovsky, 2002c): 


Rey (u, Lo, o) = R? (u, Ho, p) 
x exp[—-y(1 — cy)u(u, uo, P)]. (4.32) 


where c=0.05. The value of c was obtained by the 
parameterization of calculations with exact radiative 
transfer code, developed by Mishchenko et al., 1999. 
Kokhanovsky (2002c) found that the error of Eq. 
(4.32) is smaller than 5% at the nadir observation, 
the solar angle smaller than 75° and y < 1.18. 
Detailed studies of accuracy of Eq. (4.30) were per- 
formed by Kokhanovsky et al. (1998b). In particular, 
they found that the error of Eq. (4.30) is smaller than 
15% at t= 10 and œ> 0.95. Kokhanovsky and Roza- 
nov (2003) proposed to substitute texp(— x) in Eq. 
(4.30) by f= texp(— x) — ł, where 7 has been introduced 
in Section 4.1. The accuracy of the modified equation is 
given in Fig. 16a—d (Kokhanovsky and Rozanov, 
2003). It is better than 6% at solar angles smaller than 
70° at nadir observation, if t > 5, A < 1.5 um for the 
Deirmendjian’s Cloud C.1 model with agp=6 um. 


Eq. (4.30) can be used to determine the optical 
thickness and effective radius from spectral reflection 
function measurements over extended cloud fields 
(Kokhanovsky and Zege, 1996; Kokhanovsky et al., 
2002). 

Surface reflection is accounted for by Eq. (4.12). 
Substitution of Eq. (4.30) into Eq. (4.12) yields: 


R(u, Mo, P, T) = Roo (H, Ho, P, T) 


Irs 
exp( x y) l= rF T (u, Ho), 
(4.33) 
where 
T(u, Mo) = tKo(u)Ko(uo) (4.34) 


is the transmittance function of a cloud, r is the total 
reflectance of a cloud, ¢ is given by Eq. (4.31), and R% 
is given by Eq. (4.32). It is assumed that 
tK(W)K(uo) ~ tKo(u)Ko(uo). This assumption is valid 
for single-scattering albedos close to one. 


4.3. Total reflectance 


Let us find the approximate solution for the total 
reflectance r in Eq. (4.33). Clearly, the value of r# 
1—t. This is due to light absorption in a cloudy 
medium. 

The total reflectance or the spherical albedo r is 
defined by (Sobolev, 1972): 


1 2n T ; 
=z ao f dûsin2ð 


x f dvosin2VoR(Vo, V, P, T). (4.35) 
0 


Let us also introduce the plane albedo 


1 27 T 
R(vo) =z ao | dûsin2ðR(Vo, V, p, 1) 


and the diffuse transmittance 


1 2n T 
T(V) = =f ao f dvsin20T (Jo, V, ,T). 
T Jo 0 
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Fig. 16. (a) Dependence of the reflection function of a plane-parallel homogeneous cloud layer with water droplets on the solar zenith angle at 
nadir observation for different values of the cloud optical thickness at the effective radius of droplets equal to 6 um, the wavelength equal to 1.55 
um (exact data - symbols, approximation - lines). The droplet-size distribution is given by Eq. (2.3) at u=6 (Kokhanovsky and Rozanov, 2003). 
(b) The same as in Fig. 16a but for various effective radii and t=5, 100 (Kokhanovsky and Rozanov, 2003). (c) The error of approximation 
obtained using data in Fig. 16a (Kokhanovsky and Rozanov, 2003). (d) The error of approximation for various values of y=0.12 (1), 0.27 (2), 
0.37 (3), 0.83 (4), and 1.2 (5) and t=5 (Kokhanovsky and Rozanov, 2003). 


For the case of idealized semiinfinite nonabsorbing 
clouds (Sobolev, 1972) and as a result of the con- 
servation of energy law, 


1 


2n / 
— Vosin2v 
an Jy ao f dVosin2 vo 


n/2 
x | ddsin20R°, (Vo, V, p) = 1 (4.36) 
0 


and 


2r n/2 
— |] dọ | ddsin2VR°,(00,9,~) =1, (4.37) 
2n 0 0 


i.e., all photons injected into a cloudy medium are 
reflected back in outer space after an infinite travel 
time. Here R? (00,09,90) is the reflection function of a 
semiinfinite nonabsorbing cloud. The reflection func- 
tion R? (o0, o) of a cloudy medium only weakly 


A. Kokhanovsky / Earth-Science Reviews 64 (2004) 189-241 219 


depends on its microstructure (see Fig. 13); by defi- 
nition, it does not depend on either the optical thick- 
ness, T=Gex:L, or the single-scattering albedo, wp = 
OscalOext. Here ex, is the extinction coefficient and 
Osea 1S the scattering coefficient of a cloudy layer of 
the geometrical thickness L. 

It follows from Eqs. (4.30) and (4.35) for absorb- 
ing clouds: 


r= Tf — texp(—x—y), (4.38) 
where 

1 2n n/2 
fo => ap f dvsin2v 

2T Jo 0 

n/2 
x f dVosin2VoR... (Vo, V, P) (4.39) 
0 


and we take the normalization condition (Van de 
Hulst, 1980) 


2 f amiKo) =1 (4.40) 


into account. The approximate formula (4.10) for 
Ko(u) obeys the exact integral relation (4.40). The 
constant r „ represents the total reflectance of a semi- 
infinite layer. According to the definition (Eq. (4.36)), 
r~ =l at œ= 1. Eq. (4.39) is not readily analytically 
integrated at arbitrary values of œo, but it follows from 
Eqs. (4.18) and (4.39) as œ@ọ— 1 (see also the dis- 
cussion in the previous section) that 


To =1—y. (4.41) 
At larger values of y using the same substitution as 
was used in the derivation of Eq. (4.28) from Eq. 
(4.25), the integral (Eq. (4.39)) gives approximately 


Poo = exp(—y) (4.42) 


or (see Eq. (4.28)) r.,=/'. Similar result for the 
plane albedo reads (Kokhanovsky, 2001a): 


N. (Ho) = exp[—YKo(H0)]- (4.43) 


It easy to show that the diffuse transmittance is given 
by the following equation (Kokhanovsky, 200 1a): 


T (Uo) = tKo( uo). (4.44) 


The accuracy of Eq. (4.42) is studied in Fig. 17a,b. We 
conclude that the accuracy of this formula is better 
than 5% at wo>0.93 for crystalline clouds. For water 
clouds, the accuracy of this formula is better than 5% 
at Wo>0.97. 

Combining Eqs. (4.31), (4.38) and (4.42), we have 
for the total reflectance of a cloud layer: 


_ sinh(x + %y)exp(—y) — sinh(y)exp(—x — y) 
= sinh(x + ay) 


(4.45) 


Substitution of Eq. (4.45) into Eq. (4.33) enables 
the reflection function of a cloud over underlying 
surface with the spherical albedo r, to be calcu- 
lated. Such calculations were performed and pre- 
sented in a graphical form by Kokhanovsky et al. 
(1998b). 

Total light absorptance in a cloud layer is given by 
a=1—r-—t, where r is calculated by Eq. (4.42) and t 
by Eq. (4.31). As a result, we have the following 
analytical equation for the total light absorbtion inside 
a plane-parallel cloud having a finite thickness: 


p Sinh( + op)(1 — exp(—y)] — sinh(y)[1 — exp(—x— y) 
sinh(x + ay) 
(4.46) 

It follows as «= 1: 

—_ sinh(x) ___ sinh(y) 
sinh(x +y)? — sinh(x +y)’ 

__ sinh(x + y) — sinh(x) — sinh(y) 
a sinh(x + y) , aa) 


which yields the well-known formulae presented else- 
where (see, e.g., Zege et al., 1991). 

Overall, the global radiative characteristics of 
cloudy media are well described by only two param- 
eters, which are called the generalized radiative trans- 
fer parameters (Zege et al., 1991): 


1- Mo 
a(t): 
(4.48) 


x=7/3(1—q@)(l—g), y=4 


The parameter x describes the attenuation of a light 
field in deep layers of semiinfinite weakly absorbing 
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Fig. 17. (a) Cloud spherical albedo calculated with the exact radiative theory (Yanovitskij, 1997) (symbols) and according to approximation 
(4.42) (lines) at g=0.75, 0.8 and various wo. In exact calculations, Heney—Greenstein phase function was assumed. (b) Errors of Eq. (4.42) as 


obtained from data in Fig. 17a. 


media. For the light intensity in deep layers, [(1,t)= 
W(Wexp(— x), where the angular distribution of light 
field y(u) does not depend on the optical depth, t. We 
see from Eq. (4.42) that y=In(1/r.,). Thus, the radi- 


ative characteristics of optically thick cloud layers are 
determined by parameters x and y, which govern 
asymptotic regime and light reflection for semiinfinite 
cloudy media. 
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Table 5 

Radiative characteristics of cloudy media 

Radiative Equation 

characteristic 

Reflection R( to.) = RCh ho p) — texp(— x —y) 
function Ko(W)Ko(Ho) 


Plane albedo 
Spherical albedo 


RCo) =R CHo) — texp(— x — y)Ko(Ho) 
r=r,, —texp(—x—y) 


Transmission T( HoH) = tKo( lo) Ko() 
function 

Diffuse T(o) = tK (Mo) 
transmittance 

Global t=sinhy/sinh(ay +x) 
transmittance 


Here RB .(Ho,L.9) = RS (Host P)expl — YuCLonHP)], UoH p)=Ko 
(Ho) Ko(W)/R%(HosHs9), R o(uo) = expl[—yKo(uo)], Ko(to)= 
(3/7) (1+2uo), r..=exp(—y), #=1.07, and R°, (uou) is the 
reflection function of a semiinfinite_nonabsorbing cloud medium, 


1—00 


x= 1y/3(1 — œ)(1 — 8), y = 43T 


This concludes the consideration of reflection, 
transmission, and absorption characteristics of cloudy 
media. For the sake of convenience, we present main 
equations of this section in Table 5. They can be used 
for a rapid estimation of radiative characteristics of 
water clouds, giving that their microphysical or local 
optical parameters are known. Also they give a con- 
venient formulation and starting point for a semi- 
analytical solution of the inverse problem, namely, 
for finding t and der from spectral measurements of 
reflected light (Kokhanovsky et al., 2003). 

For instance, measurements of light fluxes above 
and below clouds were performed using pyranometers 
by Asano et al. (2002). They used look-up table 
approach to retrieve cloud properties. Taking into 
account that the cloud system studied was optically 
thick, more simple solutions given here can be equally 
applied, leading in the decrease in the computational 
burden. 


5. Polarization characteristics 
5.1. Introduction 


Cloudy media do not only absorb and scatter ra- 
diation. They also can be considered as major sources 
(along with aerosols, gases, and underlying surfaces) of 
polarized light in the terrestrial atmosphere. Clouds 
introduce a new property to incoming unpolarized 


solar light, namely, the preferential direction of oscil- 
lation of the electric vector in a light beam (Kokha- 
novsky, 2003b). Generally speaking, both reflected and 
transmitted light is elliptically polarized. 

The polarization arises mainly due to single scat- 
tering of light by a water droplet, an aerosol particle or 
an ice crystal. Molecular scattering is another source 
of the polarized light in the terrestrial atmosphere. 
Multiple scattering, which is of particular importance 
for clouds, leads to randomization both photons 
directions and polarization states. This decreases light 
polarization due to increase in the optical thickness of 
a cloud. Another effect of multiple scattering is the 
introduction of ellipticity in scattered light beams. It 
should be pointed out that the solar light scattered by 
spherical water droplets becomes partially linearly 
polarized. This is due to the fact that the ellipticity 
of singly scattered light is equal to zero for spherical 
scatterers if the incident light is unpolarized (Van de 
Hulst, 1957). 

Studies of polarization characteristics of solar light 
transmitted and reflected by cloudy media have a long 
history (Hovenier, 1971; Hansen, 1971a,b; Hansen 
and Hovenier, 1974; Gehrels, 1974; Coffeen, 1979; 
Coulson, 1988; Knibbe, 1997; Kokhanovsky, 2003b). 
However, the real burst of research in this area was 
given by a launch of the Polarization and Direction- 
ality of Earth Reflectances (POLDER) instrument on 
board the Japanese ADEOS-I satellite (Deschamps et 
al., 1994; Chepfer et al., 2001). The POLDER was 
able to transmit to ground stations a huge amount of 
information about polarization characteristics of light 
reflected from cloudy media, aerosols, and underlying 
surfaces at several wavelengths. Specifically, the first 
three components of the Stokes vector S,(/, Q, U,V) 
(Chandrasekhar, 1950) have been measured for wave- 
lengths, 2, equal to 443, 670, and 865 nm. There is no 
doubt that even more advanced polarimeters (see, 
e.g., Cairns et al., 1999) with more wide spectral 
coverage will appear on board different satellites in 
future, which makes further theoretical studies of 
polarization characteristics of cloudy media extremely 
important. This is due to potential possibilities for 
global retrievals of cloud microphysical parameters 
(Breon and Goloub, 1998), the shape of particles 
(Deschamps et al., 1994) and the optical thickness 
of clouds (Kokhanovsky, 2000) based on polarization 
measurements. 
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Of course, similar information can be obtained 
from the reflected intensity measurements. However, 
it could well appear that the degree of polarization 
can be used as a source of additional information 
about cloud particle-size distributions close to the top 
of a clouds (Breon and Goloub, 1998). Indeed, the 
high proportion of photons scattered from a thin 
upper layer of a cloud in creating light polarization 
is quite understandable. Multiply scattered light 
fluxes from deep layers are hardly polarized at all. 
Radiative characteristics, on the other hand, represent 
the cloud as a whole. Thus, the effective radius 
derived from radiative measurements is an average 
of large ensembles, possibly very different particle- 
size distributions, present in different parts of cloudy 
media. 

It has been emphasized by Hansen (197 1a) that the 
fact that the polarization data refers to cloud tops may 
be not a drawback, but an advantage. Indeed, it is 
easier to interpret observations which do not refer to 
some average values over the entire cloud depth. 

The polarization characteristics of cloudy media 
can be studied applying numerical codes (Hansen, 
1971a; Liou, 1992) based on the vector radiative 
transfer equation solution. However, one can also use 
the fact that cloud fields are optically thick in most 
cases (Kondratyev and Binenko, 1984). This allows 
application of vector asymptotic analytical relations 
(Domke, 1978) derived for optically thick disperse 
media with arbitrary phase functions and absorption. 
These solutions help clarify physical mechanisms 
and main features behind the polarization change 
due to the increase of the size of droplets or the 
thickness of a cloud. Analytical solutions also pro- 
vide an important tool for the simplification of 
the inverse problem. They can be used, e.g., in 
studies of the information content of polarimetric 
measurements. 

The drawback of well-known vector asymptotic 
theory for optically thick layers (Domke, 1978) is 
related to the fact that the main equations of this 
theory depend on different auxiliary functions and 
parameters. Calculations of these functions and 
parameters are rather complex due to their gene- 
ral applicability in terms of the single-scattering 
albedo. 

The special case of weakly absorbing light scatter- 
ing media was considered by Kokhanovsky (200 1a) in 


the framework of the modified asymptotic theory. 
Formulae obtained are rather simple and allow for 
semianalytical solution of the inverse problem 
(Kokhanovsky and Weichert, 2002). They can be 
applied to water clouds due to weak absorption of 
light by water droplets in a broad spectral range (at 
least up to 2.2 um, see Fig. 2). 


5.2. Reflection and transmission matrices 


Let us introduce the reflection R(ju,J%,~) and 
transmission Î(u,uo,p) matrices of cloudy media. 
Reflection and transmission matrices depend on the 
optical thickness, t, of a cloud, microstructure of a 
cloud, cosine uo of the solar angle Vo, cosine u of 
the viewing angle ð, and the relative azimuth y= 
W -— Wo, where Wo and y are azimuths of the Sun 
and observer, respectively. The importance of these 
matrices is that they can be used for the calcu- 
lation of the Stokes vector of reflected S,(/,Q, U, 
V) and transmitted SUI, Q’,U',V') light under 
arbitrary illumination of a cloudy layer. Namely, 
it follows for plane-parallel cloudy layers (Hove- 
nier, 1971): 


= 1 1 2n A 
S+(u, @) = F duy i do'k(u, ue — 9’) 


x So(ut, o”), (5.1) 
= a ae 
Sı(u, p) -1f duu! d do'T(u, u, p — 9’) 


— — $. 
EE E ( =) 
0 
(5.2) 


where Soli p) represents the Stokes vector of the 
incident light beam. The second term in Eq. (5.2) 
is small for most of cloudy media (tœ 1) and can 
be neglected. 

In a good approximation, the Sun is a source of 
unpolarized light incident from only one direction (Uo, 
Po). Then it follows: 


=> 


Slup) = 5 - 1) 5(W — vo) F, (5.3) 


where vector F is normalized so that the first element 
of nF (namely, nF) is the net flux of solar beam per 
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unit area of a cloud layer. Eqs. (5.1) and (5.2) simplify 
if one accounts for Eq. (5.3): 


Si(u, p) = Ru, Mo, Y — Wo) F, (5.4) 
S(t, 9) = TH, Ho, Y — Wo) F, (5.5) 


where only first element of the column vector 


Fy 


differs from zero for unpolarized solar radiation. Thus, 
only elements Rj, Tj), j= 1, 2, 3, 4) are relevant to 
studies of solar light propagation, scattering and 
polarization in terrestrial atmosphere. Other elements 
are of interest only in the case of polarized (e.g., laser) 
incident beams. 

Knowledge of matrices R and 7 is of primary 
importance for understanding radiative and polariza- 
tion characteristics of light fluxes, which emerge from 
a cloud body. They play a role of well-known Mueller 
matrices (Shurcliff, 1962; Tynes et al., 2001) for 
multiple-scattering cloudy media. Note that Eq. (5.4) 
can be written in scalar form for unpolarized illumi- 
nation: 


IT=RuFi, Q= Rafi, U= RF, V = Rai 


(5.7) 
and similarly to the transmitted light. Thus, the 
intensity of the reflected light is determined by the 
element R,, of the reflection matrix. This is the only 


element of major concern of most cloud optics studies 
up to date. The degree of polarization is given by 


P=,/P + P2, (5.8) 


where 


2 2 
A= bE (5.9) 


Ri 


is the degree of linear polarization and 
P, == (5.10) 


is the degree of circular polarization. One can show 
(Van de Hulst, 1980) that for particular viewing geo- 
metries (e.g., dg =1, “= 1, p=0,7 ) the values of R31 
and R4; are equal to zero and light is partially linearly 
polarized. The degree of circular polarization, P., is 
equal to zero in this case. 

It is common to use also a different definition of 
the degree of polarization, namely, 


Q 
P=-= 5.11 
: (5.11) 
for linearly polarized light beams. It follows from Eq. 
(5.11) at the nadir viewing geometry (R3;=R4, =0) 


that 
P=-—. (5.12) 


We will use this definition of the degree of polar- 
ization in this review. The second component of the 
Stokes vector, Q, is proportional to the difference 
I,—1,, where values of J, and J, are intensities of light 
polarized parallel and perpendicular to a meridional 
plane, which contains the normal to a layer and a 
viewing direction. Thus, light polarized perpendicular 
to the meridional plane is characterized by the positive 
polarization. Note that the sign of polarization degree 
strongly depends on the observation geometry. For 
instance, it is positive in the main rainbow region (Van 
de Hulst, 1957). It could be both positive and negative 
in the glory-scattering region. 

The degree of polarization of singly scattered 
light for solar illumination is given in Fig. 18a,b. 
They were obtained using Mie theory for polydis- 
persions of water droplets at the wavelength, 7, equal 
to 443 nm. The particle-size distribution was given 
by Eq. (2.3). 

Polarized light scattering by hexagonal crystals 
were studied by Labonnote et al. (2001) and Liou 
(2002), the first paper being devoted to a special case 
of inhomogeneous crystals. 

It should be pointed out that the Stokes vector 
elements Z, Q, U, V (see Eq. (5.7)) can be used to 
determine the polarization ellipse parameters such as 
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Fig. 18. (a) Degree of polarization of light singly scattered by a polydispersion of spherical water droplets characterized by the gamma particle- 
size distribution (Eq. (2.3)) with the half-width parameter u=6 and the effective radius aef=4, 6, 8, 10, and 16 um as the function of the 
scattering angle at 1 =443 nm. (b) Same as in Fig. la, but at a.¢=6 um and different values of the half-width parameter = 2, 4, 6, 8, 10, and 12. 


the ellipticity and the angle of preferential oscilla- Thus, reflection and transmission matrices 
tions of the electric vector in a light beam (Shurcliff, allow determination of both intensity and polar- 
1962). ization states of reflected and transmitted light 
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beams under arbitrary illumination of a cloud 
layer. 


5.3. Modified asymptotic theory 


5.3.1. General equations 

The optical thickness of water clouds is usually 
large. Some mixed and crystalline clouds are also 
extended in the vertical direction for many hundred 
meters and could be optically thick (Uchiyama et al., 
1999). Thus, asymptotic vector theory (Domke, 1978; 
Van de Hulst, 1980) is applicable to cloudy media. 
Using this theory, the azimuthally averaged reflection 
matrix Ř(u,4o,1) and transmission matrix Î(u, 40,1) of 
a homogeneous plane-parallel weakly absorbing tur- 
bid layer of a large optical thickness, t, can be written 
in the following form as y—0 (Kokhanovsky, 
2001a): 


R(t, Host) = Roo (4, Ho) — T(u, uo, t)exp(—x — y), 
(5.13) 


F (11, uo, 1) = tKo( 1) Ka (Ho), (5.14) 


where x and y are given by Eq. (4.46), ¢ is given by 
Eq. (4.31) and 


> 3 a i eee 
Rot =F | | dec? der? (ç, u, ole 
T Jo 0 
(5.15) 
where 
1 
z= . (5.16) 
0 


R.(,Ho) is the azimuthally averaged reflection 
matrix for a semiinfinite medium with the same local 
optical properties as a finite medium under study. The 
vector K” (uo) in Eq, (5.15) is transpose to the vector 
K(o). The vector Ko(u) also occurs in the so-called 
Milne problem (Minin, 1988; Wauben, 1992). It 
describes the angular dependence and the polarization 
characteristics of light emerging from semuinfinite 
nonabsorbing turbid layers with sources of radiation 
located deep inside the medium. 


Eqs. (5.13) and (5.14), derived by Kokhanovsky 
(2001a) for small probabilities of photon absorption, 
are much simpler than initial asymptotic formulae 
obtained by Domke (1978). They allow calculation 
of the reflection and the transmission matrices of thick 
weakly absorbing disperse media by simple means if 
the solution of the problem for a semiinfinite medium 
with the same phase matrix as for a finite layer under 
consideration is available (Wauben, 1992; Van de 
Hulst, 1980). This reduction of a problem to the case 
of a semiinfinite nonabsorbing medium is of impor- 
tance for radiative transfer theory, in general, and to 
cloud optics, in particular. 


5.3.2. Nonabsorbing media 

It follows from Eqs. (5.12)—(5.14) for the degree 
of polarization at the nadir observation in the case of 
nonabsorbing clouds (Kokhanovsky, 2001a) that: 


0 4 
P(u) = aa a l (5.17) 
where 
u(t) = oly (5.18) 


and ¢ is given by Eq. (4.7). Our calculations show that 
we have roughly u = 1 and, therefore (see Eq. 5.17), 
P(t) = P® (uo)/r, where we used r+t=1 for non- 
absorbing clouds. It means that brighter clouds have 
smaller values of the degree of polarization and vice 
versa. The degree of polarization for a semiinfinite 
nonabsorbing cloud at the nadir observation 


Real, Lo) 


— al 
R? i1(1, Mo) 6 ?) 


P? (Ho) = 


depends on the solar angle and the effective radius of 
particles. Here Ro iC ua) and R sitlay) are corre- 
sponding elements of the reflection matrix (see Eq. 
(5.4)) of a nonabsorbing semiinfinite medium. Note 
that Eq. (5.17) was derived accounting for the equality 
(Wauben, 1992) Ko2(1)=0. Kokhanovsky and Wei- 
chert (2002) proposed the parameterization of this 
function in terms of der at the rainbow-scattering 
region. 

Thus, Eqs. (5.17) reduces the calculation of the 
degree of polarization of light reflected by a cloud in 
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the visible to the case of an idealized situation of a 
nonabsorbing semiinfinite media. Calculations with 
Eq. (5.17) and exact vector radiative transfer equation 
are given in Fig. 19 at the rainbow-scattering geom- 
etry. The value of P? (uo) in Eq. (5.17) was obtained 
from the exact vector radiative transfer calculations. 
We see that accuracy of Eq. (5.17) is high. Note also 
the decrease of the degree of polarization with cloud 
optical thickness. This can be used for the cloud 
optical thickness determination at the rainbow-scatter- 
ing geometry (Kokhanovsky, 2000). 

In more general case (see Eqs. (5.13) and (5.14)), 
the functions Ko() and R®, (11,49,~) for a nonabsorb- 
ing semiinfinite light-scattering medium should be 
found numerically. They depend on the phase matrix 
of a random medium in question. However, by defi- 
nition they do not depend on the single-scattering 
albedo, wo, and the optical thickness, t. Another 
interesting feature of these functions is their weak 
dependence on the cloud microphysical parameters 
(size, shape, and chemical composition of particles) of 
a medium under study in the broad range of angular 
parameters. This is due to the randomization of both 


polarization states and propagation directions of pho- 
tons in highly scattering nonabsorbing semiinfinite 
layers irrespective of local optical properties of the 
medium where they propagate. 


5.3.3. Weakly absorbing media 

Eq. (5.17) is easily generalized to account for weak 
absorption of light in a cloudy medium. Namely, it 
follows from Eqs. (5.12) and (5.13) (Kokhanovsky, 
2001b): 


P( to) _ Po (uo) , 


~ Tarup r) PN 


where the global transmittance, ¢, is given by Eq. 
(4.31) and 


Ko(1)Ko(Ho) 
Ro ( 1 El Ho) 
Values of P (Ho) and R(1,u4o) represent the degree 


of polarization and reflection function of a semiinfin- 
ite weakly absorbing cloud at the nadir viewing 


u* (Ug) = (5.21) 
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Fig. 19. Dependence of the degree of polarization on the solar angle for cloud optical thickness 8, 15, and 30 at A=0.65 um, V9 =37°, dep=6 
um, = 6 (see Eq. (2.3)) according to the exact radiative transfer calculations (symbols) and Eq. (5.17) (lines). 
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geometry. They can expressed via the correspondent 
functions for nonabsorbing media. In particular, Eq. 
(4.32) can be used to find R (1,40). Similar equation 
holds for P,,(uo) (Kokhanovsky, 2002c): 


P. (po) = P®, (uo)explyu(1 — 0.05y)). (5.22) 


Eqs. (4, 32), (5.20)—(5.22) reduce the calculation for a 
weakly absorbing finite layer to a more simple case of 
a semiinfinite nonabsorbing medium. The accuracy of 
these equations has been studied by Kokhanovsky 
(2001a, 2002c). 

Note that the dependence of the degree of polar- 
ization of reflected light on the size of droplets has 
been studied by Hansen (1971a) and Kokhanovsky 
(2003b). 


6. Image transfer 
6.1. Introduction 


Let us consider now pecularities of image transfer 
in cloudy media. From vision theory, a cloud is a 
high-frequency filter. The transmission of high-fre- 
quency signals (in time and space domains) is reduced 
considerably in comparison with low-frequency sig- 
nals as a result of light scattering. Any diffused source 
of light can be considered as a superposion of point 
light sources. Thus, in linear optical systems the 
image of such an object with intensity Jo(7’) is a 
linear superposition of images of point sources (Zege 
et al., 1991): 


=f 3 dF” Io(F)S(F'.F dr", (6.1) 


where the point spread function (PSF) S(r’,7’) 
describes the process of the transformation of the 
object intensity Jo(7”) in the initial plane to the image 
intensity /(r’) in the image plane. The point spread 
function is a main notion of the image transfer theory. 
Eq. (6.1) has a more simple form in a frequency 
domain: 


I0) = SM’), (6.2) 


where 

I) = / / IF )e dF", (6.3) 
z / J DE ar”, (6.4) 

S@) = L J S(F) eo? ap", (6.5) 


and bv’ is the space frequency. One of the main 
problems of the image transfer theory is to determine 
the value of the Fourier transform of the PSF, namely, 
the optical transfer function (OTF) S(%). After that the 
value of Z (v) can be determined from Eq. (6.2) and for 
the light intensity in the image plane, it follows from 
Eq. (6.3): 


I(r") P= [i Jei” "dD. (6.6) 


The optical transfer function depends on the prop- 
erties of media between an object and an image 
(Wells, 1961; Volnistova and Drofa, 1986; Zege et 
al., 1991). Different methods of calculation of this 
function and the modulation transfer function (MTF) 


1st) 
MO) = Ts) 


are described by Ishimaru (1978) and Zege et al. 
(1991). 


6.2. Optical transfer function 


We will consider here the only case of the image 
transfer through a cloud with large particles at normal 
incidence of a light beam. In this case, the OTF is a 
real function. It does not depend on the azimuth for 
randomly oriented particles. It was obtained (Zege et 
al., 1991) that 


S(v*) = exp[—t + wotn(vx)], (6.7) 


where v* = vL is the dimensionless frequency, L is the 
geometrical thickness of a layer, v = |0| and 


noe) = f P(v*v)dy, 


P(v*v) = 7 | ` p(0)Jo(Ov*v)0d0. (6.8) 
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Here Jo(@v*v) is the Bessel function. The optical 
transfer function can be obtained from Eqs. (6.7) 
and (6.8) if one has information on the optical thick- 
ness of a cloud (< 7, in the approximations (6.7) and 
(6.8)), the single-scattering albedo, wo, and the phase 
function p(0). 

For media with the complex microstructure (e.g., 
crystal clouds), the crux of the difficulty is that the 
phase function is not known a priori. In this case, 
different models to describe this function can be used 
(see Table 6). For instance, if one supposes that the 
phase function is described by the Gaussian function: 


aie Zep |- Ua (6.9) 


where b is constant, which is inversely proportional to 
the effective radius of droplets, and @ is the scattering 
angle, it follows that 


ry)? 
P(v*v) = exp - on | (6.10) 
and, therefore, 
_ jot JT bv* 
so) = o| THz e(“)]. (6.11) 


where erf(bv*//2) is the probability integral, given 
by 


erf (=) = =| eTde 
V2) V2nJ-x , 


One can obtain at large values of frequency, v*, erf 
(bv*/\/2) = 1 and, therefore, 


(6.12) 


S(v*) = exp(—t). (6.13) 
Thus, the OTF is determined by the nonscattered 
light alone at v*œ1. It follows directly from Eqs. 
(6.7) and (6.8) as well. 

Note that at large optical thickness, the number of 
nonscattered photons is low. Thus, the value of the 
OTF is small and the contrast of the image is low. It 
follows at small values of the frequency, v*, that 


aae 


(6.14) 


Table 6 

Phase functions p(0) and their Fourier— Bessel transforms p(o) ( Y is 
the normalization constant, x is the size parameter) (Zege et al., 
1991; Kokhanovsky, 2001a) 


pO) p(o) 

(2Yexp(— Y0))/0 -t 
VI +o 

2Y7exp(— Y0) TAL +P” 


2/YPexp exp(— (Y?o7)/2) 
(— PRY’) 
oC oC oC 2 
(4s? (0x) z {aost -4V1 -8 } = 
i 
0, o > 2x 


and, therefore, 


si") = ew] f a(i eN 


+2 
M(v*) = exp - ow | : (6.15) 
In particular, we have at v* =0: 
S(0) = exp[—t(1 — @o)], (6.16) 


which means that S(0) differs from one only due to the 
light absorption process. Note that the value of b is 
inversely proportional to the average size of particles. 
Thus, one can see that the loss of contrast is higher for 
thicker clouds with larger single-scattering albedos 
and smaller droplets or crystals. 

The dependence of the OTF on the size and 
refractive index of droplets was studied in detail by 
Zege and Kokhanovsky (1994), Zuev et al. (1997), and 
Kokhanovsky (2001a) using the geometrical optical 
approximation for the phase function p(@). This al- 
lowed them to obtain analytical relationships between 
values of OTF and cloud microphysical parameters. 
Zege and Kokhanovsky (1992) used these relation- 
ships for the optical particle sizing problem solution. 

Comparisons of calculations the OTF with Eq. 
(6.7) and Monte-Carlo calculations (Drofa and Usa- 
chev, 1980) were performed by Zege and Kokhanov- 
sky (1994). It was found that the accuracy of 
approximate Eq. (6.7) for cloudy media with an 
effective radius of water droplets 6 um and t=1 is 
better than 5% in the visible. The accuracy decreases 
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with optical thickness, t, which should be less than 
5-7 to apply Eq. (6.7). 


7. Satellite remote sensing of cloudy media 
7.1. Optical thickness 


Equations presented in previous sections can be 
used for a rapid estimations of the radiative, polar- 
ization and image characteristics of cloudy media. 
They also can be used to check the accuracy of new 
algorithms using the fact that the numerical solution of 
the radiative transfer equation (Thomas and Stamnes, 
1999) and results presented above for optically thick 
layers should converge as t — œ and f — 0. 

Kokhanovsky and Macke (1999) used approxima- 
tions, presented in Table 5, to study the influence of 
the shape of particles on the radiative transfer in 
clouds. They found that clouds with nonspherical 
particles are more reflective (larger values of reflec- 
tance) as compared to clouds with spherical droplets 
with the same value of volume to particle surface area 
ratio. The opposite is true for transmittance. Similar 
results were obtained by Takano and Liou (1989). 

The most important area of the application of 
asymptotical solutions lies in the field of remote sens- 
ing (Rozenberg, 1967; Rozenberg et al., 1978; King, 
1981, 1987; Kokhanovsky and Zege, 1996; Kokha- 
novsky, 2000, 2001a; Kokhanovsky et al., 2003). 

In particular, this approach avoids or reduces (if 
thin clouds are also under consideration) the need for 
precalculation and storage of the so-called look-up 
tables, which is a standard technique in passive cloud 
remote sensing (Arking and Childs, 1985; Rossow et 
al., 1989; Nakajima and King, 1990). 

We have for the global transmittance from Eq. 
(4.14) after simple algebraic calculations: 


(l—r,)A 


= A 


where the function A is introduced. It is given by 


R? (u, uo, P) — R(H, Mo, P, 7) 
Ko(u)Ko( uo) 


A=A(u, Mo, ) = 
(7.2) 


The analytical results for functions R2(L, 110.) and 
Ko(u) have been presented above. Thus, the global 
transmittance, t, and correspondingly the total reflec- 
tance, r=1-— t, can be obtained from Eqs. (7.1) and 
(7.2) and knowledge of the surface albedo r, and the 
measured value R({1,L19,(,7). 

For such a retrieval, it is not necessary to know the 
optical thickness of clouds and the average size of 
droplets. This is an extremely important point for 
climate studies, where the global and temporally 
averaged value of the cloud reflectance, r=1 — t, is 
an important parameter. Usually, r<0.8 for natural 
water clouds in the visible (Danielson et al., 1969), 
which implies that clouds with optical thicknesses 
larger than 70—100 appear not very often (Trishchenko 
et al., 2001). The reduced reflectance in visible can be 
also related to aerosol absorption in clouds (Melnikova 
and Mikhailov, 1994, 2000) and to the inhomogeneity 
and finite size of clouds (Stephens and Tsay, 1990). 

Using Eq. (7.1), we obtain t= A at r,=0, and t=0 
and r=/ at r,=1. This shows that all photons incident 
on optically thick nonabsorbing clouds over surfaces 
with r,=1 survive and return back to outer space. 
They yield no information about actual cloud thick- 
ness. This explains why the retrieval of cloud param- 
eters over bright surfaces (e.g., snow and ice) hardly 
can be performed in visible (Platnick et al., 2001). The 
information on the global transmittance, t, can be used 
to find the scaled optical thickness (Rozenberg et al., 
1978; King, 1987), given by 


t* = 1(1—g). (7.3) 
It follows from Eqs. (7.1) and (7.2) that 
4 
i= 3 — a), (7.4) 


where ¢ is given by Eq. (7.1). Again, the value of t* 
can be obtained although there is no information 
about the size of droplets and the actual optical 
thickness of clouds. 

Eq. (7.4) can be used for the retrieval of t* from 
the measurement Å(u,u0,ọ,T) at a single wavelength 
(King, 1987). Eq. (7.3) is used for the derivation of 
the optical thickness, t, (see Eq. (7.3)) if the value of g 
is known (Rossow, 1989) (approximately 0.74 for ice 
clouds as previously discussed). However, for warm 
clouds, the asymmetry parameter g depends on the 
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size of droplets even for nonabsorbing channels (see 
Eq. (3.35)). Often, the dependence g(aer) is neglected 
and it is assumed that a.-=10 um for water clouds 
(Rossow and Schiffer, 1999). Then it follows from Eq. 
(3.35) that g=0.86 at 24=0.65 um and agp=10 um. 
Utilizing Eqs. (7.3) and (7.4), this value of g can be 
used for a crude estimation of the optical thickness of 
liquid clouds. 

Errors can be introduced if one assumes the fixed a 
priori defined value of g. It follows from Eq. (3.35) at 
2=0.65 um that g = 0.84-0.87 at age=4—20 um. From 
Eq. (7.3), we have t= S1*, X =(1—g) | = 6.3-7.6 
and t=[9.4,11.5] at t* = 1.5, depending on the value of 
g used. The assumption that aep= 10 um yields g=0.86 
and 3 =7.2, t=10.7. 

This leads to a relative error of 7—14% in the 
retrieved optical thickness (i.e., a range of possible 
values from t=9.4 to t=11.5 instead of t= 10.7). 
This uncertainty in the optical thickness can be 
removed if measurements in the near-infrared region 
of the electromagnetic spectrum are performed, ena- 
bling the size of droplets and, therefore, the asymme- 
try parameter g to be estimated. For this, however, we 
should be sure that we have a liquid and not ice or a 
mixed-phase cloud. Another uncertainty arises due to 
the possible contamination of clouds by absorbing 
aerosols (Asano et al., 2001, 2002). Then Eq. (7.1) is 
not valid. 


7.2. Size of droplets 


As was specified above for correct estimation of 
optical thickness of clouds from space, we need to 
know the effective radius of droplets. The size of 
droplets can be found if the reflection function in 
near-infrared region spectrum is measured simultane- 
ously (Nakajima and King, 1990; Kokhanovsky and 
Zege, 1996). This is due to the fact that the reflection 
function in the infrared strongly depends on the 
probability of photon absorption by droplets. This 
probability is proportional to the effective radius of 
droplets as was discussed previously (see Eq. (3.21)). 

The influence of absorption and scattering of light 
by molecules and aerosol particles on the measured 
value R(L,10,9,T) is often neglected in the cloud 
retrieval algorithms. However, correction can be easily 
taken into account if needed (Wang and King, 1997; 
Goloub et al., 2000). The influence of surface reflec- 


tion on the cloud reflection function, assuming that the 
surface is Lambertian with albedo r, is easily taken 
into account, leading to (see Eqs. (4.14) and (4.33)): 


. tı (der, w) [1 — rs1] 
Ri (der, w) = Re 
i(a fw) o0 J= ra [l — tı (der, w)] 


x Ko(u)Ko(uo), (7.5) 


Ro(der, w) = R? exp(—y(aet)(1 — cy(aer))u) 
- { exp(—n(aer») — ya) 
to (er, W)I's2 \ 


1— Trs2r2(def, w) 
X t(aef, w)Ko(u)Ko(uo)- (7.6) 


The subscripts “1” and “2” refer to wavelengths, 2, 
and A>, in visible and near-infrared channels, respec- 
tively. The values of r,ı and 7. give us the surface 
albedos in visible and near-infrared respectively. The 
explicit dependence of functions involved on the 
parameters aef and w to be retrieved is introduced in 
brackets. The liquid water path w is preferred to the 
optical thickness in retrieval procedures due to the 
independence of w on the wavelength. The optical 
thickness is uniquely defined if agp and w are known. 

Eqs. (7.5) and (7.6) form a nonlinear system of 
two algebraic equations having two unknowns (der 
and w), which can be solved by standard methods 
and programs. In particular, we can find the value of 
w from Eq. (7.5) analytically (Kokhanovsky et al., 
2003). Substitution of this result in Eq. (7.6) gives us 
a single transcendent equation for the effective radius 
of droplets determination (Kokhanovsky et al., 
2003). The accuracy of this semianalytical retrieval 
algorithm has been studied by Kokhanovsky et al. 
(2003). 


7.3. Thermodynamic state of clouds 


The discrimination of liquid water from ice clouds 
is of importance for many applications including 
flight safety and Earth climate studies. The size and 
shape of particles in warm and ice clouds are different. 
This influences the energy transmitted and reflected 
by a cloud. 

This discrimination can be performed, taking into 
account the difference in angular or spectral distribu- 


A. Kokhanovsky / Earth-Science Reviews 64 (2004) 189-241 231 


tion of reflected light. We present results of calcula- 
tion of the reflection function of cloudy media with 
liquid and frozen water droplets in Fig. 20. It follows 
from this figure that minima in the reflection function 
of ice clouds (e.g., near 1.5 and 2.0 um) are moved to 
larger wavelengths as compared to the case of liquid 
droplets. Of course, this is due to the difference in 
spectral behaviour of the imaginary parts of the 
complex refractive index of liquid water and ice (see 
Fig. 2). Minima for liquid water also moved to larger 
wavelengths as compared to the absorption bands of 
water vapour. These different positions of minima can 
be easily registered with modern spectrometers (see, 
e.g., Dvoryashin, 2002; Knap et al., 2002). 

Another possibility is to consider different angular 
behaviour of the reflection function for ice and water 
clouds at specific scattering geometries (e.g., rainbow, 
glory and halo scattering). In particular, the reflection 
function of water clouds, as distinct from ice clouds, 
has a maximum near the rainbow-scattering angle, 
which also can be easily detected. This feature 
becomes even more pronounced if the degree of 


reflection function 


“0.0 0.5 1.0 
wavelength, um 


polarization at the rainbow geometry is studied 
(Goloub et al., 2000; Kokhanovsky, 2003b). 


7.4. Cloud-top height and cloud fraction 


Another important characteristic of a cloud is its 
height. It can be retrieved using data from space-borne 
lidars (Winker and Trepte, 1998). Passive measure- 
ments also can be used. For instance, Yamomoto and 
Wark (1961) proposed the use of the oxygen A band, 
centred at 0.761 um. The physical basis of this 
method depends on a deep minimum around 0.761 
um due to oxygen absorption. This minimum is not 
shown in Fig. 20 because only scattering and absorp- 
tion of light by cloud particles was accounted for in 
the calculation for this figure. The depth of the 
absorption line will depend on the cloud height. 
Photons can hardly penetrate thick clouds and be 
absorbed by the oxygen in the air column below the 
cloud. This will increase the value of the reflection 
function at 0.761 um for the case of clouds at high 
altitudes. The depth of the absorption line is larger for 


1.5 2.0 2.5 


Fig. 20. Spectral dependence of the reflection function of cloudy media for the nadir observation and the solar angle equal to 60°. Clouds are 
composed of water or ice spherical particles with effective radius 6 um. It is assumed that particles in a cloud are characterized by the gamma 
particle-size distribution (Eq. (2.3)) with the parameter .=6. The geometrical thickness of cloud is equal to 500 m. The liquid water path is 100 
g/m”, which gives the optical thickness of 27 at the wavelength 0.55 um. The reflection of light from surface and the scattering and absorption of 


light by aerosols and gases are neglected. 
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low clouds. Practical application of the method, how- 
ever, is not so simple (Kuze and Chance, 1994; 
Koelemeijer et al., 2001). First of all, the depth of 
line also depends on the oxygen absorption cross 
section, which varies with temperature and pressure. 
Thus, one should use a priori assumptions on the 
temperature and pressure variation with height in the 
Earth’s atmosphere. The generally unknown surface 
albedo and cloud geometrical thickness also can 
influence the retrieval accuracy. Other possible sour- 
ces of errors are described in detail by Kuze and 
Chance (1994) and Koelemeijer et al. (2001). 

The largest complication arises for pixels that are 
only partially covered by clouds. Ignoring them will 
lead to a big reduction of data; hence, to overcome 
this problem, Koelemeijer et al. (2001) proposed an 
algorithm that retrieves simultaneously effective 
cloud-top height/pressure and cloud fraction, assum- 
ing that cloud is a Lambertian surface. Such an 
assumption leads to underestimation of cloud-top 
heights as compared to in situ measurements. 

Global information on cloud fraction/cover, Q, is 
of a considerable importance by itself (Batey et al., 
2000). Usually, the value of Q varies in the range 
0.55-0.85, depending on the exact region under 
study. Globally, clouds cover around 60% of the 
atmosphere. This once more underlines the impor- 
tance of clouds in the radiation balance and atmos- 
phere heating rates studies. 

Palle and Butler (2002) state that the global cloud 
cover increased during the past century. They argue 
against a dominating role of the solar activity (via 
galactic cosmic rays) on the cloud formation. 


7.5. Remote sensing of crystalline clouds 


Remote sensing of optical thicknesses and effective 
size of droplets in crystalline clouds is complicated by 
low optical thickness (usually smaller than 5), the 
high spatial and temporal variability of cloud proper- 
ties, nonspherical shape of particles and their possible 
horizontal orientation. If the thickness of a crystalline 
cloud is high, then the optical thickness can be derived 
in the same way as was discussed above, assuming an 
asymmetry parameter of 0.74 (see Eqs. (7.1)—(7.4)). 
The problem is to find the corresponding reflection 
function of a semiinfinite medium (see Eq. (7.4)). One 
possibility is to calculate it beforehand using the phase 


function of a fractal fictive scatterer. Mie theory 
should not be used in this case. 

For thinner clouds, one should build the precalcu- 
lated table of reflection functions, which should be 
compared with experimental data to establish both the 
optical thickness of clouds and the size/shape of 
droplets (Masuda et al., 2002). This is not an easy 
problem. In particular, the model of ice spheres cannot 
be used for this purpose (Mishchenko et al., 1995, 
1996; Chepfer et al., 1998; Rolland et al., 2000; 
Doutriaux-Bucher et al., 2000). 

Yang et al. (2001) considered the influence of 
vertical variability of ice crystals habits and crystal 
dimensions on the satellite ice cloud retrieval algo- 
rithms. They accounted for the fact that more complex 
shapes and larger sizes of crystals are more typical for 
the base of the cloud as compared to its top, where 
crystals are smaller and more rounded. The authors 
state that the vertical distribution of optical character- 
istics of clouds can be neglected if visible channels are 
used, but the vertical inhomogeneity should be fully 
accounted for if one is interested in the average size of 
particles retrievals (e.g., from the measurements of the 
cloud reflection function at 2.11 um moderate reso- 
lution imaging spectrometer channel (King et al., 
1992)). Information on the vertical structure of a 
crystalline cloud is not known a priori, which compli- 
cates the retrieval procedure. We see, therefore, that the 
creation of a suitable look-up table for reflection func- 
tions of crystalline clouds is not at all a trivial problem. 

Note that ice clouds are often studied using 
ground-based systems. This allows us to derive the 
cloud optical thickness and the effective radius of 
particles. For this, measurements of direct light and 
light scattered at small angles usually are performed. 
The important problem here is to account for multiple 
light scattering using the radiative transfer equation 
(Thomas and Stamnes, 1999). For layers with optical 
thickness smaller than approximately 5, a so-called 
small-angle approximation (Zege et al., 1991) can be 
used. The intensity of the diffused light, Za, is given by 
a simple integral in this case (Zege et al., 1991; 
Kokhanovsky, 200 1a): 


naD) = f O - exp(—2)Lo(od)ade, 


(7.7) 
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where (see Table 6) 
I(t) = exp[—(1 — @op(o))] (7.8) 


is the normalized mutual coherence function for 
electric fields separated by the distance o/k for a plane 
incident wave, where k=2z/A; Ip is the intensity of 
incident light and it was assumed that a cloud is 
illuminated along the normal. The function 


(a) = 5 |, POOO (7.9) 


is the Fourier—Bessel transform of the phase function 
p(0). The angle V gives the angle that measured 
transmitted light beam makes with the normal to a 
cloud layer. Most frequently used functions p(@) and 
their transforms p(o) are given in Table 6. Function 
(7.8) can be used to study the loss of coherence of 
optical waves in cloudy media (Zege and Kokhanov- 
sky, 1995). 

In conclusion, we note that both ice and water 
clouds are often studied with lidar techniques (Platt, 
1978a,b; Bissonnete et al., 1995; Sassen, 1991, 2001; 
Winker and Trepte, 1998; Mishchenko and Sassen, 
1998; Zege et al., 1998a,b; Winker, 1999; McGill et 
al., 2002). 


8. Inhomogeneous clouds 
8.1. Vertical inhomogeneity 


Vertical and horizontal inhomogeneity of clouds 
can be dealt with in the framework of the Monte- 
Carlo methods of the radiative transfer equation 
solutions (Schreier and Macke, 2001). This is of 
considerable importance for crystalline media, as was 
discussed in the previous section. Unfortunately, 
Monte-Carlo methods are extremely slow and cannot 
be applied in the operational satellite cloud retrieval 
algorithms. 

Accounting for the vertical inhomogeneity, how- 
ever, can be easily done in the framework of the 
theory of optically thick layers as discussed previ- 
ously. Corresponding equations are presented by 
Sobolev (1972) and Yanovitskij (1997). For the sake 
of simplicity, we consider here only the case of clouds 
in the visible, where light absorption can be neglected. 


Then the reflection function of a vertically inhomoge- 
neous optically thick cloud can be presented by Eq. 
(4.6) as in the case of a homogeneous layer. However, 
the meaning of parameters in this equation is different 
in an inhomogeneous case. That is, t is the optical 
thickness of a vertically inhomogeneous cloud, R Lis 
the reflection function of a semiinfinite cloud with the 
same vertical distribution of optical characteristics as 
a finite cloud, and g is the average asymmetry 
parameter. The value of g does not change consid- 
erably with droplet size (see Eq. (3.35) as Xe — 00); 
thus, it can be assumed to be a priori defined value 
(say, e.g., 0.86; Rossow and Schiffer, 1999). Then 
accounting also for the weak sensitivity of the reflec- 
tion function R° to the size of droplets (see Fig. 13), 
we state that Eq. (4.6) can be also applied to the 
derivation of optical thickness of vertically inhomo- 
geneous clouds. This suggests that one cannot retrieve 
the vertical distribution of extinction coefficient in a 
cloud from the reflection function measurements in 
the visible. Lidar returns, however, can provide this 
information (Zege et al., 1998a,b). 

The situation in the near-infrared is not so simple 
because the derived size of droplets depends on the 
average photon penetration depth, which is of course a 
function of wavelength (Platnick, 2000). It means that 
the derived effective size of droplets is a function of 
the wavelength for vertically inhomogeneous clouds. 
Generally, the absorption increases with wavelength; 
hence, penetration depth decreases with wavelength. 

Thus, we arrive to the conclusion that reflection 
functions at larger wavelength (e.g., 2.11 um) will 
give us values of the droplets radii closer to the cloud 
tops and reflection functions at smaller wavelengths 
(e.g., 1.6 um) give us the radius of droplets deep 
inside clouds. These different values of radii are 
considered as a shortcoming of the retrieval method 
for vertically inhomogeneous clouds. On the other 
hand, this opens a new possibility to study the vertical 
distribution of droplets’ sizes, analysing spectral 
reflectances of clouds, obtained both from airborne 
and satellite platforms. 

In conclusion, we note that radii of droplets 
increase with height (Brenguier et al., 2000). It means 
that values of radii retrieved at 1.6 um should be 
smaller than those obtained from the channel at 2.11 
um. This is also observed experimentally (Platnick et 
al., 2001). 
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8.2. Horizontal inhomogeneity 


Real cloud fields are horizontally inhomogeneous 
(Titov, 1998; Marshak et al., 1995, 1998) and have 
complex shapes, which complicates theoretical mod- 
elling (Rogovtsov, 1991, 1999). This also produces 
unphysical angular dependence of cloud optical 
thickness for different solar angles (Loeb and 
Davies, 1996; Loeb and Coakley, 1998) if the model 
of the horizontally inhomogeneous layer is used in 
the retrieval procedure for inhomogeneous cloud 
fields. Fouilloux et al. (2000) found that derived 
cloud optical thicknesses and effective radii for 
inhomogeneous clouds depend on the averaging 
scale. Therefore, they state that comparisons between 
aircraft measurements and satellite observations are 
not valid for heterogeneous clouds (which is the case 
for most clouds). This puts the validation of cloud 
satellite products with airborne radiometers in ques- 
tion. 

Approximations for horizontal photon transport 
within real-world clouds were developed by Platnick 
(2001). In particular, he derived analytic approxima- 
tions for the root-mean-square horizontal displace- 
ment of reflected and transmitted photons relative to 
the incident cloud-top location. 
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Usually, the influence of horizontal inhomogenity 
of clouds on their radiative characteristics is studied in 
the framework of the independent column (or pixel) 
approximation (Cahalan et al., 1994). This approxi- 
mation neglects the horizontal photon transport 
between adjacent columns. Then one can obtain for 
the reflection function of a pixel in the framework of 
this approximation (Pinkus and Klein, 2000): 


k= f R(t)f (t)dz, (8.1) 


where f(t) is the optical thickness distribution function 
for a given pixel and R(t) is the reflection function for 
a horizontally homogeneous cloud with a given opti- 
cal thickness, t. Cahalan et al. (1994), Fu et al. (2000), 
and Schreier and Macke (2001) found that this 
approximation has a high accuracy for the domain- 
averaged radiative fluxes. 

It is known (Cahalan et al., 1994) that R< R(t), 
where 7 is the average optical thickness, which means 
that the values of 7 obtained from measurements over 
horizontally inhomogeneous clouds in the assumption 
of a horizontally homogeneous plane-parallel clouds 
are underestimated. Thus, correction of the optical 
thickness obtained by an empirical factor is needed 
(Rossow and Schiffer, 1999). 
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Fig. 21. The coefficient of variance of the cloud optical thickness spatial distribution (Barker et al., 1996). 
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It is difficult to apply Eq. (8.1) to the radiative 
transfer problem analysis due to the necessity to 
perform calculations R(t) many times for a given 
pixel. The problem can be greatly simplified (Kokha- 
novsky, 2003d) if one uses approximate analytical 
solutions, which are valid for optically thick clouds 
(see e.g., Section 4), in combination with analytical 
forms for the function f(t), which is usually given by 
the gamma (Barker et al., 1996) or the lognormal 
(Nakajima et al., 1991) distributions. The values of 
the coefficient of variance of this distribution obtained 
from experimental measurements over a number of 
cloud scenes are given in Fig. 21, which was prepared 
using data of Barker et al. (1996). We see that the 
coefficient of variance is in the range 20-80%. 

The lower limit of integration in Eq. (8.1) should 
be 10 to allow for the use of Eq. (4.14) in integral 
(8.1). The integral in the range from 0 to 10 can be 
estimated using the exact radiative transfer equation. 
This contribution is small at large values of t and can 
be neglected. 


9. Summary 


The goal of this review was to consider various 
approaches to calculate local and global optical 
characteristics of cloudy media. Main tools reviewed 
are the geometrical optics approximation and the 
theory of optically thick turbid layers. Both theories 
can be substituted by exact solutions (the Mie theory 
for local optical characteristics and the radiative 
transfer theory for the radiative characteristics) in 
the case of homogeneous clouds with spherical par- 
ticles. 

However, they appear to be very important in 
bringing the forward propagation models closer to 
the reality (e.g., account for nonspherical shape of 
crystals, effects of cloud inhomogeneity). Approxima- 
tions developed allow consideration of cases that are 
difficult or impossible to handle with exact techni- 
ques. Another important feature of the approximate 
methods is the possibility to simplify inverse prob- 
lems of cloud optics. 

The cloud optics studies are far from being com- 
pleted. Most results of this review are valid only for 
homogeneous plane parallel clouds, which is a good 
starting point for cloud optical properties investiga- 


tions. In reality, most clouds have very complicated 
shapes and structures. Accounting for these features 
will be a main topic of cloud optics research for years 
to come. The progress in respect to the description of 
local optical characteristics of crystalline clouds with 
particles of complex shapes is expected due to large 
efforts currently undertaken in this direction. 

Other topics of high priority, which will be devel- 
oping in a future, are studies of clouds with lidars 
(e.g., space-based). This will allow provide informa- 
tion on cloud geometrical parameters with a high 
spatial resolution. Depolarization and polarization 
techniques will also gain further momentum in cloud 
optics studies. 

All of these developments will allow us to charac- 
terize cloud properties and their temporal and spatial 
variability on a scale that was not possible in the past. 
In turn, the information obtained will make it possible 
to understand and describe the influence of clouds on 
climate and its change, which remains a major unre- 
solved problem of a modern climatology. 
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